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TOPOLOGICAL PROPERTIES OF DIFFERENTIABLE 
MANIFOLDS* 


BY HASSLER WHITNEY 


I. INTRODUCTION 


1. The Problems. In many fields of work one is led to the con- 
sideration of n-dimensional spaces. A given dynamical system 
has a certain number of “degrees of freedom” ; thus a rigid body, 
with one point fixed, has three. A line in euclidean space is de- 
termined by four “parameters.” We therefore consider the posi- 
tions of the rigid body, or the straight lines, as forming a space 
of three, or four, dimensions.t But when we try to determine 
the points of the space by assigning to each a set of three, or 
four, numbers, we are doomed to failure. This is possible for a 
small region of either space, but not for the whole space at once. 
The best we can do is to cover the space with such regions, de- 
fine a coordinate system in each, and state how the coordinate 
systems are related in any two overlapping regions. They will 
be related in general by means of differentiable,{ or analytic, 
transformations, with non-vanishing Jacobian. Any such space 
we shall call a differentiable, or analytic, manifold. 

For a complete study of such spaces, we must know not only 
properties of euclidean n-space E”, which we may apply in each 
coordinate system separately, but also properties which arise 
from the manifold being pieced together from a number of such 
systems. It is these latter properties, essentially topological in 
character, which form the subject of the present address. 

Suppose we wish to study differential geometry in the n-di- 
mensional manifold M”. At each point p of M*, the possible 
differentials (or “tangent” vectors) form an n-dimensional vec- 


* An address delivered before the Society on September 9, 1937, in State 
College, by invitation oi the Program Committee. 

t The first space forms the group G* of rotations in 3-space; it is homeo- 
morphic with projective 3-space P*. The second space is homeomorphic with 
the total space (see §2) of the tangent vector space of the projective plane P?, 
or of the 2-sphere S* if we use oriented lines. This is easily seen by considering 
together all parallel lines. 

t Differentiable will always mean continuously differentiable. 
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tor space V(p), the so-called tangent space at p. For topological 
considerations, it is sufficient to consider vectors of unit length,* 
or directions, which form a sphere S(p) of dimension n —1. This 
set of spheres forms the tangent sphere-space of M”. Most of our 
work will be on the problem, how do the spheres fit together 
over the whole manifold? Suppose M™* is imbedded in a higher 
dimensional manifold M™ (for instance, euclidean E”). Then 
we may consider the normal unit vectors at each point, form- 
ing an (m—n—1)-sphere, and thus the normal sphere-space. The 
methods we use in Part II were discovered independently by 
E. Stiefel [11] and myself [15]. Here we shall not attempt 
always to describe the most elegant methods, but lean rather 
to the intuitive side. 

In the last part we give some fundamental results, due to 
de Rham,f{ on the theory of multiple integration in a manifold. 
If we wish to integrate over an r-dimensional subset, and no 
measure function is given, we integrate “differential forms,” in 
other words, “alternating covariant tensors” of order 7; we shall 
call these simply r-functions. It is necessary for various pur- 
poses to find what “exact” r-functions exist. This problem, in 
any small region, is rather trivial; again, we must consider the 
whole manifold to answer it. The methods we describe are some- 
what different from those of de Rham. 


2. Sphere-Spaces. A sphere-space S(K) is defined as follows. 
Let Ss be the unit v-sphere§ in E’+!. Let K be a complex. To 
each point p of K let there correspond a v-sphere S(p); if p¥q, 
we assume that S(p) and S(g) have no common points. For each 
closed cell o of K and each point # of a, let &,(p) be a (1-1)map 
of So into S(p); let &.(p, g) be the point into which the point q 
of S¢ goes. If the cells ¢ and a’ have the common point p, then 


* As every differentiable manifold may be imbedded in a euclidean space 
(see the bibliography, Whitney [14], Theorem 1), we may define in it a Rieman- 
nian metric and thus lengths of vectors, and so on. 

+ Numbers in square brackets refer to the bibliography at the end of this 


paper. 
t See [6], [7], and [8]. The main theorem was suggested by E. Cartan, 
Comptes Rendus de |’Académie des Sciences, Paris, vol. 187 (1928), pp. 
196-198. 
§ That is, the set of points )-x? =1 in E’. A 0-sphere is a pair of points, 
a 1-sphere is a circle, and a 2-sphere is a spherical surface. 
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&.(p) and &.(p) are both defined. We assume that 


which is a (1-1) map of So into itself, is orthogonal,* and varies 
continuously with p. We call K the base space, and the set 
S(K) of all points on all S(p), the total space.j If K were a gen- 
eral point set, we would replace the cells o by open sets covering 
K. We might replace the S(p) by vector spaces V(p) (in the 
obvious manner for tangent and normal spaces); however, the 
topological properties of these are reducible to those of sphere- 
spaces. For any sphere-space S(K) and subset L of K, there is 
a corresponding sphere-space S(L), the part of S(K) over L. We 
shall call the map £,(p), (pin a), the o-coordinate system in Sic). 
It determines (for a fixed orientation of So’) one of two possible 
orientations of each S(p), (p in a). 

We find invariants which serve to distinguish between sphere- 
spaces with the same base space by studying the following ques- 
tions. First, is it possible to choose, for each p in K, a point 
o(p) of S(p), so that ¢(p) is continuous in K? Any such map 
we shall call a projection of K into S(K). More generally, is it 
possible to find k projections ¢:(p), - - - , dx(p) of K into S(K), 
so that these are orthogonal in each S(p), (1SkSv+1)? If so, 
we call S(K) k-simple. If S(K) is (v+1)-simple, we call it simple. 

For tangent sphere-spaces, the existence of a projection is ob- 
viously equivalent to the existence of a continuous field of non- 
vanishing (tangent) vectors. The only closed orientable surface 
on which there is such a field is the torus.{ (If we try, on the 
sphere, making all vectors point north or east for instance, this 
fails at the two poles; see §6.) Note that if we can find one such 
field on an orientable surface, we can at once find a second field 
of orthogonal vectors; at each point, looking along the first vec- 
tor, we let the second vector point to the left. It follows that 
the tangent space is simple in this case. 


* It is sufficient to assume that the map is linear. However, this general case 
is easily reducible to that given. Because of the assumption, we may define 
orthogonality on any S(p). 

+ For a study of some total spaces, see Hotelling [4], Seifert [9] and [10], 
Stiefel [11], and Threlfall [12]. 

t For a study of this problem in manifolds, see Hopf [3]; also Alexandroff- 
Hopf, Topologie I, pp. 548-552, and Stiefel [11]. 
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Suppose S(K) is simple. Let gi=(1, 0,---, 0),--+, gous 
=(0, 0,---, 1) be the unit points of Sy’. Set ¢(p, ¢:) =4<(p), 
(i=1,---,v+1). Then for each p in K, we can define {(p, q) 
for all other g in So” uniquely so that it is an orthogonal map of 
So” into S(p). Clearly ¢ is (1-1) and continuous. Hence if S(K) 
is simple, then S(K) is the cartesian product of K and S,’. 

If the tangent space to M™ is simple, then to each direction 
¢(p, at p corresponds a direction ¢(p’, gq), (same q), at any other 
point p’. Thus a parallelism may be defined throughout M” if and 
only if the tangent sphere-space is simple. 

Suppose the normal space to M* in M™ is simple. Then corre- 
sponding to each normal vector function ¢,(p) we may define a 
function f;(p) in M™ with ¢,(p) as its gradient at p, and vanish- 
ing in 1", We may then define M” as the set of points in M™ at 
which all the f;(p) vanish. Conversely, if it is possible to define 
M™ in this manner, with the gradients of the f;(p) linearly inde- 
pendent, then the normal space is simple.* 


II. INVARIANTS OF SPHERE-SPACES 

3. Elementary Properties of Complexes. Let "07 be the inci- 

dence number of the cells ¢7 and ¢;7~'. Define the boundary and 
coboundary of the chain A’ =)_a,o; by 


r r—1 r r+1 3 r+l 
(2) 0A = 0; 6A = a0; 


Then A’ is a cycle or cocycle if 0A’ =0 or 6A’ =0; A’ is homologous 
to B’, A’~B’, if A*—B’ is a boundary; similarly, for cohomolo- 
gous, for which we use the symbol —. If we identify homologous 
r-cycles, we obtain the rth homology group; similarly for the rth 
cohomology group. Define scalar products of r-chains by 


(3) ao): Bio’) = >, aB;. 


Note that A’-o,7 is the coefficient of o7 in the chain A’. As 
and both equal the incidence number 


* See Seifert [10], and Whitney, Annals of Mathematcis, vol. 37 (1936), 
pp. 865-878. 

t For further details, see for instance Whitney, Matrices of integers and 
combinatorial topology, Duke Mathematical Journal, vol. 3 (1937), pp. 35-45, 
and On products in a complex, Annals of Mathematics, vol. 39 (1938). We shall 
refer to these papers as MI and PC respectively. 
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we have, on equating, multiplying by coefficients, and summing, 
(4) Br! = Ar-6Br, 


We shall use o also for the set of points in the closed cell ¢, and 
Oo for its point set boundary. Let K* be the r-dimensional part 
of K, the sum of all cells of dimension <r. 


4. Orientability of Sphere-Spaces. First note that the part of 
S(K) over K®, S(K®°), is simple. We can define S(K®) as a prod- 
uct by setting ((a) = £,(a) for each vertex a. We now ask, under 
what conditions is S(K') simple? We attempt to extend ¢(p) 
through the 1-cells of K. Consider any 1-cell o=ab. S(a) and 
S(b) are given orientations both by ¢ and by &,. If the orienta- 
tions given by ¢ are both the same as, or both opposite to, those 


Fic. 1 


given by é,, set F'(¢) =0; otherwise, F'(0) =1. In Fig. 1,* if the 
end point of each vector, say at a. is a’={(a, qi), we have 
F'(ab) = F'(bc) =1, F'(ca) =0. We may say, F'(c) is 0 or 1 ac- 
cording as the orientations given to S(a) and S(b) by ¢ are simi- 
lar or opposite, when viewed through the o-coordinate system. 
We define a characteristic 1-cocycle of K by letting F'(a) be the 
coefficient of ¢ in a chain: 


(5) F' = coefficients integers mod 2. 


* For p in ab, all vectors £:(p) point up or all point down, 
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To show that F' is a cocycle, we shall find the coefficient 
6F!-¢? = F!-d0? of 6F! in any 2-cell o? of K. For each 1-cell ab 
of dc”, the similarity or dissimilarity of the orientations given 
to S(a) and S(6) when viewed through the ab-coordinate system 
is clearly the same as when viewed through the o?-coordinate 
system. Say 00?=ab+bc+ --- +ea. If we compare the orienta- 
tion of S(a) with that of S(b), then of S(b) with that of S(c),..., 
the number of changes of orientation is even, as we come back 
to the original orientation of S(a). Hence F!-d07=0 (mod 2), 
as required. 

We might have chosen different coordinate systems over K', 
and have thus found a different F'. To find the change made in 
F', suppose the orientation of S(a) is changed,* while that of 
each other S(x) is unaltered. Then for each 1-cell o with a asa 
vertex, F'(a) is changed from 0 to 1 or from 1 to 0, that is, it is 
increased by 1 (mod 2). Thus the new F' is the old plus 6a 
(mod 2). Thus we may change F! only by coboundaries, and 
the cohomology class h*(F') is an tnvariant of S(K). Clearly if F! 
is a characteristic cocycle, we may obtain any other F’!— F! 
by reorienting some of the S(a) for vertices a. 

We call the cohomology class h!=h'(F‘) the characteristic 
1-class of S(K). If it vanishes, we call S(K) orientable, otherwise 
non-orientable. Suppose S(K) is orientable. Then we may re- 
orient some of the S(a) so that F!=0. For each 1-cell ab, S(a) 
and S(b) now have either both the same or both the opposite 
orientations with the a- and b-coordinate systems as with the 
ab-coordinate system. In the latter case, change the orientation 
in the ab-coordinate system so that the orientations agree. 
Carrying out the same process in the 2-cells, and so on, we 
have finally chosen all coordinate systems so that for any point 
p of K in two cells o and a’, their two coordinate systems give 
S(p) the same orientation. This justifies the term “orientable.” 

To answer our original question, we shall show that if S(K) 
is orientable, and only then, it is simple over K'. If ht=0, choose 
the £&, so that F!'=0 and set ¢(a) = £,(a). Take any 1-cell ab. By 
hypothesis, 


* That is, let 2 denote a reflection in So’, and set &’(a) = &.(a)2. Clearly 
F' is independent of the o;'-coordinate systems. 
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are orthogonal transformations of S¢ into itself, that is, points 
¥(a) and (bd) of the group G’+! of orthogonal transformations in 
(v+1)-space. As F1(ab)=0 (mod 2), ¥(a) and y(b) are in the 
same component of G’+! (of which there are two), and hence we 
may join them by an arc. Let ¥(p) run along this arc as p runs 
along ab, and set 


(7) 5(p) = 


in ab. Thus ¢ is defined in ad, and in the same manner, through- 
out K!. Conversely, if ¢ is defined in K'!, set £,(a) =¢(a) for each 
vertex a; this makes F!=0. 

Let A be any 1-cycle (mod 2) of K; then A- F! is a numerical 
invariant (mod 2) of S(K). For if we had used F’! instead of F', 
then F’!— F!= 6G° for some G®, and 


(8) A-F" — A-F! = A-8G® = 0A-G = 0 (mod 2). 


Any closed path in K defines a 1-cycle A; the orientation of 
S(K) is “preserved” or “reversed” on going around the path ac- 
cording as A - F'=0 or 1 (mod 2). 

It is easily seen that the tangent space to a manifold is orient- 
able if and only if the manifold is orientable; the orientation of 
S(K) is reversed on going around a path if and only if the orien- 
tation of the manifold is reversed. For an M” in an M™, the 
1-class of the tangent space of M” plus that of the normal space 
gives that of the part of the tangent space of M™ over M” 
(mod 2).* For an example, consider a Mébius strip in E*. 

In the future, unless the contrary is specifically stated, we 
consider only orientable sphere-spaces. 


5. Higher Dimensional Invariants. One can find, for oriented 
sphere-spaces, a characteristic cohomology class h” for 2<r 
<v+1, in the same manner as we defined orientability. We 
shall illustrate it with y=1, r=2; for instance, for the tangent 
space of S?. As S(K") is simple, we may define ¢(p, q) for pin K'. 
Set qi); this is a projection of K! into S(K!). (Fora 
general vy, we would use v orthogonal projections in K'.) 

Now take any 2-cell o of K, and consider the coordinate sys- 
tem in S(a). As pruns around = £.(p, gi) runs around 
acurve C=y(0c) on the torus S(0c) (-ee Fig. 2). Now consider 


* Compare Seifert-Threlfall, Topologie, pp. 272-276, 
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¢o(p). As p runs around @o, it runs around a curve D=¢(00), 
cutting C in the positive sense (given by the orientation of S(p)) 
say a times; it is equivalent to running around C once, and a 
times around an S(p). Set F?(¢) =a, and 


(9) 72 =: F2(¢?)o? coefficients integers. 


This chain is easily seen to be independent of the &,. (the ori- 
entation of the S(p) is preserved). Just as before, we show that 


Fic. 2 


it is a cocycle, and its cohomology class is an invariant of S(K); 
also, we may obtain any F’?— F?, 

For a general v and 7, we may always choose v —r+2 orthog- 
onal projections of K’—! into S(K*—!); we then consider these 
maps in G@(o’), in relation to £,,, to determine F’(o’) and thus F’. 

Suppose (see Stiefel [11], §6) we replace the first projection 
¢i(p) used in defining F’, for each p, by the diametrically op- 


} 
feo! 
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posite point of S(p). It turns out that if r is odd, the new F’’ is 
But F’"— F’, hence 2F"—0, and 2h” =0. We find the gen- 
eral theorem: 


For each r,2Srsv-+1, there is a characteristic r-class h’, which 
is an invariant of S(K); the coefficients used are integers unless 
r<v+1 and r is even, in which case we reduce mod 2. If r is odd, 
2h" =0; S(K) is (v—r+2)-simple over if and only if h”"=0. 


In regard to h’, see §7. 

For a non-orientable sphere-space, h’ is an invariant if we re- 
duce mod 2 in all cases.* 

As in §4, we can define a numerical invariant A’- F’ corre- 
sponding to any cycle A’. Further, if we do not reduce mod 2, 
then to any cycle A’ (mod yp), that is, such that 0A”’=yB, 
corresponds the invariant (mod yp), A’- F’; for if F’" is another 
characteristic cocycle, 


(10) At-F’" — = = = 0 (mod 
We discuss the possibility of further invariants in Part IT. 


6. Tangent and Normal Sphere-Spaces. We shall consider 
(a) the highest dimensional invariant F’*! for tangent spaces, 
(b) the same tor normal spaces, (c) relations between the spaces, 
and (d) other questions. 

(a) In a closed orientable M”, there is a “fundamental n- 
cycle” Z", whose multiples give all m-cycles. A single -cell o” 
forms a cocycle, whose multiples ko” determine all cohomology 
classes.f Any n-chain A” is ko” with k= A”-Z"=sum of coeffi- 
cients of A”. Hence, to find the characteristic n-class of S(M/*), 
we need merely find Z”- F’. 

Let us do this for the tangent space of the 2-sphere S?. The 
equator E cuts S? into the hemispheres oy and os. Orient these 
so that doy = —00s=E; then Z?=a0n+<0s. To define F?, we must 
choose a projection in K!=£, that is, a point ¢(p) of S(p) for 
each p in E. Let ¢(p) be the direction at p pointing south. To 
find F?(onv), project oy stereographically onto the tangent plane 


* Better, the theorem still holds if we replace chains in K by chains in an 
abstract complex K* obtained as follows. The cells of K* are those of K. Set 
+70, (if the latter is 0)‘ according as the o;’- and 
systems determine the same or opposite orientations of S(p) for p in o;7—. 

+ Compare MI, Appendix I, 


— 
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Py at the north pole N;* then the ¢(p) go into directions ¢’(p) 
pointing away from N. As we run around the projection of E, 
these directions make one complete turn in the same direction, 
that is, in the positive sense as given by oy. Hence F?(cy) =1. 
Similarly, projecting onto the tangent plane at S, we find 
F?(¢s) =1. Hence 


(11) Z?-F? = F*(ovy) + = 2 0. 


It follows that a continuous field of directions on S? is impossi- 
ble. 

Similarly. for the n-sphere S", we find Z"- F"=2 if 1 is even 
and =0 if m is odd.j For any closed M” (orientable or not), 
Z". F* is its Euler-Poincaré characteristic x.{ By definition of 
the n-class, this equals the “index sum of singularities” of a vec- 
tor field in M". As 2h" =0 for odd, it follows that 2x =0, x =0, 
for 2 odd, as is well known. 

If 1/" = M? is closed and orientable, then h? vanishes also, and 
it follows easily that the tangent sphere-space to any closed ori- 
entable 3-manifold is simple.§ 

(b) Consider first the normal space to a surface (say a Mébius 
strip) in E*. Let K be a subdivision of the surface. To find F', 
we choose at each vertex a one of the two normals aa’. Suppose 
that for each 1-cell ab, we join the end points a’ and b’ of the 
respective normals by an arc which lies close to ab (see Fig. 1). 
Then clearly F'(ab) is the number of intersections (mod 2) of 
a’b’ with the surface. Thus F! may be found as follows. Deform 
the 1-dimensional part K'! of K slightly so that the new positions 
a’ of the vertices a are off the surface; then F'(ab) is the 
“Kronecker index” (mod 2) of the deformed a’b’ with the sur- 
face. 

More generally, consider an M” in an M™ with m<2n; if 
m>n-+1, we assume that both are orientable, so that the nor- 


* Taking So' in Py, we may define &., as follows. Any direction dat a point 
p of cx projects into a direction d’ in Py, and determines thus a point d* of 
So'; we set Eey(p, d*) =d. 

+ It is not known whether or not the tangent space of S” is simple for all 
odd 2. 

t See (for the orientable case) Stiefel [11], §5, or Alexandroff-Hopf, 
Topologie 1, p. 549. The associated complex K* of M” has a fundamental 
n-cycle even if M* is non-orientable. 

§ See Stiefel [11], §5. 


| | 


1937-] DIFFERENTIABLE MANIFOLDS 795 


mal space is orientable. To find a characteristic cocycle F”~", 
take a subdivision K of M”, and deform the (m —n)-dimensional 
part into a new position K’"-" in M™ so that K’™~*—! does not 
touch Then is the Kronecker index (0”~"- 
of o”™-" with M” in M™”; the coefficients are integers unless 
m=n-+1. For any chain A™-", is the Kronecker 
index (A”~"- M"), as this is true if A”~" is an (m—n)-cell. 

As an application, suppose M” is a closed orientable manifold 
in E”. By joining each point of M” to a fixed point in E”, we 
construct a “singular chain” Y*+! bounded by M” (or more 
properly, by the fundamental n-cycle Z" of M"). Set s=m—n, 
and 


(12) Gu = 1)" 


Then, by a simple property of Kronecker indices (see for in- 
stance Lefschetz, Topology, page 169, (20)), 


=G*"- = (— 1)*>> (of-1- 
i 


(- 1)*(d02 -Y"+") = (of = F*-o,* 


and it follows that F"-"=6G"-"-!, and h™-"=0. 

Consider two particular cases. In m =2n, the normal space is 
1-simple, and it follows that there is a continuous field of nor- 
mals to any closed orientable M” in E?".* If m=n+2, then the 
2-class vanishes; as the S(p) are 1-dimensional and the normal 
space is orientable, it follows easily that S(/”) is simple. 

We remark that any h’ may be studied for M* in E”, 
(m>n-+r), by projecting into E”+’.{ Further. using methods 
in [14], we may construct a closed orientable M‘ in E® contain- 
ing an S? with S?- F?=1 (mod 2); thus lower dimensional in- 
variants may be non-vanishing. 

(c) For an S? in E* or in E‘, the tangent space has the in- 
variant 2, while the normal space is simple. With this in mind, 
we shall look only for relations between the invariants of the 
two spaces after reducing mod 2. A complete relation for orien- 


* The theorem is probably not true in the generality stated in [14], foot- 
note 25. 

+ See Seifert [10], Satz 1, and Whitney [15], §8, (d). 

t See Whitney [15], §9. 
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tability was given at the end of §4. For an orientable M* in an 
orientable M™, it may be shown that the same statement holds 
for the characteristic 2-classes, after reducing mod 2. This is 
probably not true for the non-orientable case. (Consider a pro- 
jective plane in E*.) For still higher dimensional invariants, no 
complete results are known. But if the normal or tangent space 
is simple, the duality theorem holds for any dimension.* It fol- 
lows immediately, by results in (a) and (b), that a closed orient- 
able manifold M” with odd characteristic x cannot be imbedded 
in E"**.+ This is the case, for instance, for the complex projec- 
tive plane, for which x =3. 

(d) Consider two sphere-spaces S*(K) and S’(K) with the 
same base space K. There is a unique sphere-space S“t’*!(K), 
their product, in which the first two spaces may be imbedded so 
that for each p in K, S*(p) and S’(p) are orthogonal great 
spheres in S*t+’+!(p). As an example, for 1" in M™, the product 
of the tangent and normal spaces gives the part S"~!(M™) of the 
tangent space of M™ over M”. A sphere-space which may be 
expressed as a product is reducible. A simple sphere-space is com- 
pletely reducible, that is, it is a product of (v+1) 0-sphere-spaces; 
any even number of these may be taken as non-orientable, the 
rest being simple. A k-simple sphere-space is the product of k 
simple 0-sphere-spaces and another space. See also §9. 

The normal space to a closed orientable M” in Et! is of 
course simple. But, one may ask, in what manner does this 
space lie in the tangent space of E**!? The outward normal at 
each point p of M” determines a point ¥(p) of the unit sphere 
So" in E**+!; thus M* is mapped by y into So’. The degree of y 
is the generalized curvatura integra of M” in E”*'.f If 1 is even, 
it is independent of the imbedding of 1/” in E”*'!, and equals 
one half the characteristic of M”. 

For another problem on the imbedding of one sphere-space 
in another, consider a manifold JJ” with a boundary B in a 
manifold M4”. Suppose a vector distribution is given over B; 
can it be extended over /"?§ We may also consider sets of vec- 


* See Stiefel [11], §6, No. 2. 

+ See Seifert [10], Satz 2. 

t See H. Hopf [2], for various results in the direction noted. 
§ See Morse [5]. 
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tors, and so on. Of course, in defining the invariants, we consid- 
ered exactly this problem, with M* replaced by a single cell oc. 
The general problem of discussing how one sphere-space lies in 
another is almost untouched. 


III. GENERAL PROPERTIES OF SPHERE-SPACES 


7. On the Classification of Sphere-Spaces. We shall now study 
in how far the invariants characterize a sphere-space with a 
given base space. But first we must clear up a matter of defini- 
tion. We have allowed changes of coordinate systems in a 
sphere-space, considering the space to be unchanged thereby. 
We shall say that two sphere-spaces S(K) and S’(K) with the 
same base space K are equivalent if there is a continuous function 
f(b) in K, such that, for each p, f(p) is an orthogonal map of 
S(p) into S’(p). In other words, there is a homeomorphism be- 
tween 6(K) and 6’(K) in which S(p) and S’(p) correspond for 
each p, the homeomorphism being orthogonal there. Then if we 
identify corresponding poinis of corresponding spheres, one 
sphere-space is obtained from the other merely by changing 
coordinate systems. If S(K) and S’(K) are both oriented, and 
there is an f(p) which preserves orientations, we call them posi- 
tively equivalent. 

If K is an open or closed cell, then its cohomology groups 
vanish, so there can be no invariants as above described. But 
also any S(K) is necessarily simple. It will turn out that the in- 
variants characterize sphere-spaces also whenever v =0 or 1, and 
whenever dim (K) $3. 

If v=0, there is just one type of S(K) for each 1-cohomology 
class of K with coefficients mod 2. Suppose now that v=1 and 
S(K) is oriented; if not, we use the methods in the footnote to 
§5. To show that S(K) is characterized by its 2-class, suppose 
that S’(K) has the same 2-class; we must prove S(K) and S’(K) 
to be positively equivalent, that is, we must find the map f of 
S(K) into S’(K), preserving the orientation of the S(p). We 
first choose projections $;(p) and (p) of K' into S(K') and 
S’(K') so as to determine the same F? = F’*. As S(K) is oriented, 
we may choose a projection ¢2(p) of K' into S(K') orthogonal 
to the first, so that, for each p, the ordered pair ¢:(p), $2(p) de- 
termines the positive orientation of S(p). Do the same for S’(p). 
These projections determine coordinate systems {(p) and £’(p) 
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throughout K! with ¢(p, g;) =$.(p), (¢=1, 2), and similarly for 
set 


Now take any 2-cell o. For each p in da, set 


(14) Wp) = EF 


this maps S¢ into S(p), then into S’(p), and then back into S/, 
and is a rotation. Thus ¥(p) is a point of the group G? of rota- 
tions of S, into itself; topologically, G? is a circle. It is easily 
seen that F’?(¢) = F?(a) implies that the degree of y is 0.* Hence 
we may extend y throughout o so that it is continuous, and de- 
fine f there by inverting (14): 


(15) = (p)- 


Thus f is extended throughout K*. The boundary of a 3-cell is 
a 2-sphere, and any map of a 2-sphere into a circle G? may be 
shrunk to a point;t hence the same process extends f through- 
out K?, and so on. 

Next take any v, and suppose dim (K) $3. The above proof 
applies again, with G? replaced by G’t!. That f may be extended 
through each a? follows from the fact that any map of a 2-sphere 
into any G’+! may be shrunk to a point.{ What has happened 
to the 3-class, which should differentiate further between 
sphere-spaces? The above proof shows that it can play no role. 
In fact, given the 2-class, the sphere-space (over K*) is com- 
pletely determined, and hence so is the 3-class.§ 

In the next simplest case, y= 2 and dim (K) =4, the invariants 
are insufficient. To show this, note that for K = S‘, the cohomol- 


* Define C and D in S(@0), using £, and ¢, as in §5, and define C’ and D’ 
similarly. Under f, D goes into D’, while C goes into C* say. As F-() = F’(o), 
C cuts D, and hence C* cuts D’, the same algebraic number of times that C’ 
cuts D’. Hence C* cuts C’ zero times algebraically, ¥(p. qi) cuts q: in So' zero 
times algebraically, and the degree of y is 0. 

+ See for instance Alexandroff-Hopf, Topologie 1, p. 516, Satz II. 

t This does not seem to have been published yet. It has been proved inde- 
pendently by E. Cartan, Ch. Ehresmann, and myself. 

§ The exact relation between h? and h' is as follows. As F? is a cocycle 
(mod 2), if we choose G? with integer coefficients so that reducing mod 2 
gives F?, then 0G?=2F* for some F*; then F* is a characteristic 3-cocycle. 
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ogy groups of dimensions 1, 2, and 3 all vanish, so that there 
can be no invariants of the above sort. Now let o; and a2 be the 
4-cells into which S* is divided by a great 3-sphere r. To define 
abstractly an S(.S*), we suppose that coordinate systems £,, and 
£,, are given, and we shall define the relation between them, that 
is, the orthogonal maps $(p) =£7,'(p)é.,(p) of Sé into itself 
(pin 7). We must thus choose a map ¢ of the 3-sphere 7 into the 
group G*. As the latter is homeomorphic with projective 3- 
space, whose covering space is the 3-sphere, these maps may be 
chosen in an infinity of distinct ways, thus defining an infinity 
of non-equivalent sphere-spaces. To define a 4-dimensional in- 
variant distinguishing between them, we must use, for instance, 
the “homotopy groups” of maps of a 3-sphere into a 2-sphere 
as coefficient group in the 4-chains considered. 

To classify sphere-spaces with vy = 2 and dim (K) =4, more ad- 
vanced methods must be used.* 


8. An Imbedding Theorem. We shall show that all possible 
sphere-spaces can be generated from certain simple ones. If 
S(K) is a sphere-space and K’ is any complex (or point set), 
then each continuous map ¢ of K’ into K generates a sphere- 
space S’(K’) as follows. For each p’ in K’, let S’(p’) be the 
sphere S(¢(p’)). 

The spaces S|v, u| of great v-spheres in the y-sphere So} form a 
set of universal sphere-spaces; any S’(K) may be generated by 
mapping K suitably into the base space of S[v, uw], with 
uw=v+dim (K).f We shall illustrate the proof for y=1. For each 
vertex a of K, let ¢(a) be an arbitrary point of the base space of 
S[1i, uw]. Let f(a) be an orthogonal map of S(a) into the corre- 


* Such methods have been found, but the classification has not been car- 
ried through. Similar methods were used in classifying the maps of a 3-complex 
into a 2-sphere; see Whitney, this Bulletin, vol. 42 (1936), p. 338. 

+ Note that the space S[1, 2] of (non-oriented) great circles on S,? is 
equivalent to the tangent space of P?; the points on each circle give the direc- 
tions from one of the corresponding poles. The base spaces of the S|», »] 
have been studied recently by C. Ehresmann, Journal de Mathématiques Pures 
et Appliquées, vol. 16 (1937), pp. 69-100. If we consider only oriented sphere- 
spaces, we may use the space of oriented great v-spheres in So. 

t If K is a compact metric space of finite dimension, we may easily prove 
the theorem by imbedding K in some E?, defining a sphere-space with an open 
set in E? containing the image of K for base space, and applying the theorem 
to this space. 
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sponding sphere S(#(a)). We next extend ¢ and f over any 1-cell 
ab of K. If we define f(p, qi) and f(p, ge) for p in ab so that they 
are orthogonal points of So” for each p, and attach to the given 
values at a and J, then there is a unique orthogonal map f() of 
S(p) into a great circle of So“ which takes on the values already 
given; this defines ¢(p) also. 

We may let f(, g:) run along any curve on So“ from f(a, gi) to 
f(b, q:) as p runs along ab. Now for each # in abd, there is a great. 
(u—1)-sphere S*-\(p) of points on orthogonal to f(p, qi); 
(2, @) must lie on this sphere. These spheres form a sphere- 
space S*~!(ab), which is simple, as ab is a 1-cell (see §7). With 
the help of a coordinate system ¢{(p, q) in it, it is easy to join 
to f(p, gz) by the required arc, as 121. The extension 
of f through cells of higher dimension is carried out in the same 
manner. 


9. On Simple Sphere-Spaces. lf S(K) is simple, it can be ex- 
pressed as a product K X Sy’, that is, a coordinate system ¢(, q) 
may be defined for all in K. In how many ways may we do 
this, preserving the orientation of the S(p)? To answer this, let 
Co be a fixed coordinate system. Then any other one, ¢, deter- 
mines a rotation $(p) ={—'(p)fo(p) of So’; there is thus defined 
a continuous map ¢ of K into G’t!. Conversely, any ¢ deter- 
mines a new coordinate system ¢ from the given one {o. Thus 
the expressions of S(K) as a product, preserving orientation, corre- 
spond to the maps of K into the group G’*' of rotations in (v+1)- 
space. 

For an example, take v=1, K a circle; S(K) is a torus. G? 
is a circle, and there are an infinity of non-homotopic maps and 
corresponding coordinate systems. For y=1 and any K, the es- 
sentially distinct expressions of S(K) as a product correspond to 
the elements of the 1-dimensional cohomology group of K with 
integer coefficients. 


IV. r-FUNCTIONS IN A DIFFERENTIABLE MANIFOLD 


10. An Imbedding Problem; Vector Functions. We begin by 
considering the problem of imbedding a differentiable mani- 
fold M* in euclidean space. Suppose f is a (continuously) dif- 
ferentiable map of M” into E”. Then f(p) has m coordinates 
filb), --- .fm(p), and if we use a coordinate system (x, , Xn) 
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in a region of each xn) is differenti- 
able. The ith contravariant coordinate vector at p is 
mapped into a contravariant vector in E™ with components 
Of:/Ox:,---, Ofm/dx;. If these m vectors, (¢=1,---, m), are 
independent, then clearly a neighborhood of » goes into E™ in 
a (1-1) way. If this holds for each p in M", we call the map 
regular. It is regular at p if the matrix 


ae 
a ith row = map of ith coordinate vector, 
Ox, Ox, 

of, Of m ; 

jth column = gradient V/; of f;, 

OX, 


is of rank n. The jth gradient Vf;, a covariant vector in M", 
may be considered as the “dual map” into M™* of the jth co- 
variant coordinate vector of E™. A regular map 1s always possible 
if m=2n (see [14], Theorem 3). To study the question for 
m<2n, perhaps the best method is to look for m covariant 
vector functions in M” such that (a) they are gradients, and 
(b) at each point of M™, n of them are independent.* 

We are thus led to the problem, when is a given (covariant) 
vector function v a gradient? As is well known, v is a gradient 
if and only if its integral {,v=f,)_v.dx; around every closed 
curve C in M* vanishes. Now let us ask, when is a vector func- 
tion such that integrating around any curve gives the same re- 
sult as integrating around any sufficiently nearby curve? The 
answer is, if and only if the integral around any sufficiently 
small curve is 0; or, the function is exact, that is, all number: 
;; in any coordinate system in M* vanish. 
In any simply connected region, a vector function is a gradient 
if and only if it is exact; but this does not hold in general over 
the whole manifold. For the relation of these functions to chains 
in a subdivision of the manifold, and generalization to higher 
dimensions, see the next section. 

For a thorough study of gradients Vf and their singularities 


* If we use [14], Theorem 4, we need not distinguish between contravari- 
ant and covariant vectors. 
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(critical points of f), we refer to the work of Morse and others.* 
Almost no study has been made yet of sets of gradients. 


11. r-Functions in a Manifold.t By an r-function in a differ- 
entiable manifold M” we shall mean a differentiable§ alternating 
covariant tensor @ of order r. We may find the integral fa of a 
over an r-cell ¢; then we define the integral over an r-chain by 


(16) a= a. 
of 

Corresponding to any r-function a there is a derived (r+1)-func- 
tion a; a is exact if a’ =0. A 0-function is a scalar function; as its 
derived is its gradient, it is exact if and only if it is a constant 
(at least if M” is connected). For a vector function, that is, a 
1-function, with components v;, the derived has the components 
W,;; as given above. The function a is derived if it is the derived 
of an (r—1)-function B. As (a’)’=0, any derived function is 
exact. By Stokes’s theorem, for any r-function a and any (r+1)- 
chain A, 


for this holds for any a’. This corresponds to (4). 
Take any subdivision K of M*. To each r-function a corre- 
sponds an r-chain, with real numbers for coefficients, defined by 


(18) o(a; f a, o(a) = o(a; 
then {,a=A -¢(a). By (7) and (4), for any 
o(a’)-o = fe a = $(a)-do = 
da 


and hence, as an equivalent of Stokes’s theorem, 


(19) o(a’) = d(a). 


* See for instance M. Morse, The Calculus of Variations, Ch. VI. 

+ See A. E. Pitcher, this Bulletin, vol. 43 (1937), p. 167. 

t Compare de Rham [6], [7], and [8], and, for non-continuous functions 
with general coefficient groups, J. W. Alexander, Proceedings of the National 
Academy of Sciences, vol. 21 (1935), pp. 511-512. 

§ The theory holds if we use tensors with continuous partial derivatives of 
order k, that is, of class C*, the derived tensors being of class C*. 
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From this formula it is apparent that if a is exact, then o(a) is 
a cocycle, and if B=a’' is derived, then $(8) is a coboundary. 

We shall need an existence lemma: If a is an exact r-function, 
and B is an (r—1)-chain such that 6B =(a), then there is an 
(r—1)-function B such that ¢(B) =B and B’ =a. We define B over 
K cell by cell, using an extended form of a lemma of de Rham* 
and a theorem on the extension of differentiable functions.t 
From this we deduce two facts: (a) If a is exact and ¢(a) is a 
coboundary, then a is derived. For B exists. (b) For any cocycle 
B there exists an exact function 6B with ¢(6)=B. For we may 
use a=0. 

We now prove the fundamental theorem: Jf we identify any 
two exact r-functions whose difference 1s derived, we obtain the rth 
cohomology group of K with real coefficients. We map each class 
of exact functions a into the cohomology class of ¢(a@); this is 
easily seen to be a well defined homeomorphism. That it is ac- 
tually an isomorphism follows at once from (a) and (b) above. 

Each cocycle A may be considered as a linear function L 
defined on cycles, with LB) =0, if we set L(B)=A-B. The 
group described is isomorphic with the set of all such functions. 
Hence, if Ci,---, C, form a maximal set of independent 
cycles (p is the rth Betti number of the manifold), the integrals 
J¢, may be chosen arbitrarily, a being determined up to de- 
rived functions. These integrals are called the periods of a. 


12. On Products of Functions.t Given an r-function a and an 
s-function 8, we define as usual their (alternating) product, giv- 
ing an (r+s)-function a®. We recall that 


(20) (aB)’ = a’B + (— 1)’af’. 


Also, chains of dimensions r and s may be multiplied, giving a 
chain of dimension r+s; if certain simple conditions are satis- 
fied, including 


(21) 5(Au B) = 6AUB + (— 1)"Au SB, 


* See de Rham [6], §21, Lemma II; we do not need his Lemma III in our 
method of proof. 

+ Whitney, Transactions of this Society, vol. 36 (1934), pp. 63-89, Lemma 
2 

t Compare de Rham [6], §27, and [7]. The proof below clearly holds for 
both orientable and non-orientable manifolds. 
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then there is a uniquely defined corresponding product in the 
cohomology groups.* We wish to show that these products cor- 
respond. 

First, it is possible to constructf an r-function w(o’) corre- 
sponding to each r-cell o” (for all 7), such that w(o") vanishes 
outside the star of and setting w(> =) aw(o;7), we 
have 


(22) w(A) = (6A),  ¢(0(A)) =A, (all A). 


Also, if J is the sum of the vertices of K, then w(J) =1. 
Now to each r-chain A and s-chain B we construct a corre- 
sponding (r+s)-chain 


(23) Au B = $[w(A)w(B)]. 


Applying (19) and (20), we derive (21). Also 7u B=B. Noting 
also a certain local condition, we see that it follows that the 
product thus defined has the required properties. We now prove 
easily in turn: (a) If @ is exact, then so is w(¢(a)), and 
$[a—w(¢(a))]=0, so that a—w(¢(a)) is derived. (Use (a) of 
$11.) (b) If B is exact also, then aB—w(¢(a))8, and also 
aB—w(d(a))w((B)), are derived. (Use (20).) Finally, 


(24) $(aB) = o(a) u 


Hence the products of functions and the products of cocycles de- 
termine the same products in the cohomology groups. 
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ON THE REMAINDER IN THE APPROXIMATE 
EVALUATION OF THE PROBABILITY IN 
THE SYMMETRICAL CASE OF JAMES 
BERNOULLI’S THEOREM* 


BY C. D. OLDS 


1. Introduction. In this paper we consider the symmetrical 
case of James Bernoulli’s theorem in the theory of probability. 
We let m represent the number of successes of an event in a 
series of m independent trials with constant probability p= 
1—gq=1/2 for the success of each trial. Then we seek the prob- 
ability P of the inequality 

| n | 

(1) 

where € is a given arbitrary positive number. The probability P 
is usually given by an approximate formula without mention of 
the error term or remainder involved.f In 1926, D. Mirimanofff 
discussed this error term and gave results which are similar, but 
not as free from restrictions as those obtained here by entirely 
different methods. § 


2. The Exact Expression for P. Let T,, represent the probabil- 
ity of m successes in the m trials and consider its generating 
function 


>» m 


m=0 


where ¢ is an arbitrary variable. It has been shown|| that 


* Presented to the Society, April 3, 1937. 

t See, for example, I. Todhunter, A History of the Mathematical Theory of 
Probability, 1865, pp. 548-553. 

t D. Mirimanoff, Le jeu de pile ou face et les formules de Laplace et de J. 
Eggenberger, Commentarii Mathematici Helvetici, vol. 2 (1926), pp. 133-168. 

§ The author wishes to acknowledge the assistance rendered him by Pro- 
fessor J. V. Uspensky. 

|| For this and similar results see A. A. Markoff, Wahrscheinlichketts- 
rechnung, 1912, pp. 18-44. 
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m=0 


In this last expression we set t=e**, multiply through by e-‘¢, 
and then integrate between the limits —7z and 7; we find 


1 
Tn = 


since 


0, msn, 
= 
m=n. 


Now let «= —}+(n/4)"¢ and express the inequality (1) in the 
form 


where /, and /, are integers. Then the probability P has the exact 
expression 


l, l, 1 
m=l, m=l, 
1 n l, 
= — (cos *) 
2 
Using the known identity 
i 


m=l, 
2 


and substituting the values of /, and /2, we find that 


1 5 n sj 1 1/2 
(2) P= (cos 
0 


sin — 
2 


3. Three Lemmas. Let \ be an arbitrary number such that 
<7. We use the expansion* 


* L. L. Smail, Elements of the Theory of Infinite Processes, 1923, p. 245. 
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3 


2 1 2 
log cos x = (2 1) r Bix? + - (2 1) ~ Box 


- 


where B,, Bo, B3,--- are the Bernoullian numbers. Conse- 
quently all the coefficients in this expansion are positive. Hence, 
we can deduce 


2 
3 — s— =—.+ 
(3) og cos > ; + Mo 

¢* 
4) — log cos — = — + — N¢q‘, 
(4 og cos + No 
where 


sec — — =4, 
2 8 


O0<NS (log -—- =) = 5, 
2 8 192 


provided 
Likewise, from the expansion* 


— 2) k 


2 sin — 
? 


where again all the coefficients are positive, we find that 


(5) =— +— 
sin — 
where 
2 
2 12 


ided On 


*K. etalk Theory and A pplication of Infinite Series, 1928, 
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Using (3) we can easily show that 


n 
(cos *) — 
2 


From (4) we find that 


(6) < nagie O< 


(cos (no? (— nNo* + 
2 192 


where 0 <£<1, and consequently 


| } ng? 
@) 2 192 


4. Application of (5), (6), and (7) to (2). Applying (5) to the 
integral on the right of (2), we have 


0 2 


1 
(cos *) “p-sin 


127 0 


1 n 
+ —f (cos -L-§-sin 
T 0 2 


1 n sj 
+—f (cos *) 
us 


sin — 


For brevity, we shall let the integrals in (8) be denoted by J;, I2, 
I3, and I, respectively. 
The inequality (7) shows that 


h= ~) dg + Aq), 
192 ¢ 


where 


2 10246 + 12288a? 
0 


an? 
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The integral in (9) splits into two integrals, the first of which 
gives 


0 
(10) in 
SIN 
us 0 
where 
2 8 
Aw | < (ne e n/ 
since 


due? 
f —< s > 0. 
u x? 


The integral in the right member of (10) is 


0 A 


us v 


et 
(=) f e~ [2dy 
0 


We replace the second integral from (9) by 


e (mo 1893 sin + Acs), 
where 
Away | < = i (1 ~) (nd2)/8 
Jy 3an 8 
We have 


as is clear if we differentiate three times with respect to a the 
integral 


__| 
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1 1/2 
(11) e~*"* cos (Bu)du = a>Q0, 
0 a 


and make obvious substitutions. 


5. The Integral I;. lf we apply the inequality (6) to J, we get 


1 
I, = e (no) sin + Avs), 
where 
Also, 
1 
e~ sin do 
127 0 
1 
where 


1 1 
| Avs) | < — = —— 
127 


3an 
Using (11) again we find that 
127 J 


6. The Integrals I; and I,. For the integral J; we have 


0 2 0 


Likewise for J; we have 


sin — 
y) 


since, for 


[December, 


Also, 
which shows that 
4 
I4| < /8 
ny? 


7. The Remainder A. Conclusion. Combining the above re- 
sults we find that the probability P of the inequality 


1 
4 
is given by 


0 


where for the remainder or error term A we have 


8 4 1 w 
< | (2+ don) | 


nh? 3nr 
with 


128 
o= Ze a + 12288a? + 10246 + 32c. 


We now let \ take on numerical values between 1 and 2, and 
seek, corresponding to each, the smallest value of m such that 
|A| <(n-*), 2. It is found that, by taking \=1.8, we have the 
inequality 


1 
nN 


for n=17. 


STANFORD UNIVERSITY 
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THE SUCCESSIVE ITERATES OF THE STIELTJES 
KERNEL EXPRESSED IN TERMS OF THE 
ELEMENTARY FUNCTIONS* 


BY D. V. WIDDER 


1. Introduction. The Stieltjes kernel is the function 


1 
Ay (é, n) 
We define its successive iterates by the recurrence relation 
H,- n dt 
0 


That these integrals all exist will appear from later considera- 
tions. Simple computation shows that 


dt log — log n 
-{ > 0,7 > 0). 


It is natural to inquire if it is also possible to express the higher 
iterates in terms of the elementary functions. It is the purpose 
of the present note to prove that this is the case. We show, in 
fact, that H,(£, 7) is a linear combination of the functions 


(log — log n)?**! (log  — log n)?* 
(2k + — n) (2k) + 7) 
the constants of combination being the coefficients of the power 


series expansion of (7/sin ms)". The precise result to be proved 
is contained in Theorem 2 of this paper, stated as follows: 


(k = 0,1,2,---), 


THEOREM 2. If 0<E< x, 0<n<~™, then 
= A 2n,2k [log é _ log n|?4-! 


0) = ; = 1,2,---), 
[log — log ]?* 
Honsilé, = 
on+1(&, 7) (2)! 


(n = 0,1,2,---), 


* Presented to the Society, September 10, 1937. 


| 
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where the constants A,,,»-x are defined by the expansion 


ms 
(1) ( ) = + + + 


sin 


2. Iterates of the Stieltjes Kernel. We introduce the successive 
convolutions of the function sech(t/2) as follows: 


hi(x) = sech (x/2), 


(2) h,(x) = f — t)dt, (n = 2,3,---). 


Since h(x) is an even function of class L, it follows that h,(x) 
has the same property. The existence of the successive integrals 
is obvious. 

We are able to express the function H7,,(£, 7) in terms of h,(x) 
as follows: 


(3) H,,(e7, e%) = (24) tw) (x — y). 


It is a simple matter of verification to check this result for n =1. 
To prove it by induction, first assume it true for »—1. Then 


fH, We, 
H,,(e7, e”) -f du 


a + e* 
= u 
(2x) ("2/2 
= - sech —— h,_1(u — y)du 
Qnelzty)/2 


—2 < 


= (27) (*-?) f hy(t)htn_i(x — y — 


= (2x) [2}—np— (2+ (x y). 


This completes the induction. 


3. Several Lemmas. In order to establish our principal result 
we first state several lemmas. 


LemMA 1. For all real x 


f e'**h,(t)dt f sech(t/2)dt 


—20 


= (27)!/? sech rx. 


| 


1937-] ITERATES OF THE STIELTJES KERNEL 
This is a familiar result.* 
LEMMA 2. For all real x 
(2) f (t)dt = sech 3x)” 


This follows from the known facts concerning the Fourier 
transform of the convolution of several functions. 


3. If 0<x<1, m=1, 2, 3,--- , then 
m! 


1-—e* (x + 
The integral on the left converges absolutely if 0 <x <1 since 
{™ 
= oe), (I> «), 
= Ole"), (I+ -«), 
for every positive e. Then 
f ———~ dt -{ = > 
t = m t 
i—et 
m! 
(x R) m+1 


k=—20 


k=1 0 


The term by term integration is justified by dominated conver- 


gence. 
Lemma 4. If —«<t<o and m=1, 2, 3,---, then 
1 1/2+ix m!' 
= e*! —_—— ds. 
dri 1/2—ize k=—20 (s + 


This follows from Lemma 3 by a general inversion formula 


for the Laplace transform. f 
* See, for example, R. E. A. C. Paley and N. Wiener, Fcurier Transfcrms in 


the Complex Domain, 1934, p. 41. 
t See, for example, S. Bochner, Vorlesungen iiber Fouriersche Integrale. 


Leipzig, 1932, p. 148. 
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LemMA 5. If 0<x<1, m=0,1,2,--- , then 


— (— 1)*mn! 1)*m! 
f = lim 

If m is positive the proof is similar to that of Lemma 3 and is 
omitted. If m=0, 


ez N 
f lim (— dt 
J 0 N-= 


N 
tlim >> (— 


Np-Nt 


ll 
M= 
| 
| 
| 


Since 
| — — (— 91 + 

< te), OSt< a), 
and since the right-hand side of this inequality is integrable in 


the interval (0, «), the interchange of limit and integration 
signs is justified. 


6. If <t<x, m=0, 1, 2,---, then 


im 1 al 1)*m! 


The proof is the same as that of Lemma 5. 


4. Evaluation of the Successive Iterates. We shall first evaluate 
the functions h,(x) in terms of the elementary functions. The 
‘result is contained in the following theorem: 


THEOREM 1. For all real x the functions h,(x) defined by (2) 
have the expression 


2 (2n—1)/2 x 
(4) he, (x) = (=) csch — > A on, 2k 
us 


(2k — 1)! 
4, 


| 
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(5) = (=) sech > 
2 imo (2k)! 
(n = 0,1,2,---). 
By use of Lemma 2 and by the known reciprocal formulas for 
the Fourier transform* we have 


h,(x) = (on) f at)"dt, (- <x < 


By the change of variable s = (1/2) —it this becomes 


(6) = f e7*(csc ms)"ds, 
1/2—iso 
(-x 
Now if the constants A,,,-; are defined by (1), it is clear that 
= 0, (k = 1, 1,72; +). 


The Mittag-Leffler theorem gives us 


A nn A n,n—-1 
(3 csc 75) (- 1) +b + (s 


k=—20 
(7) 
Aaa 
+ 
stk 
with the understanding that 
A N A 
— = lim 
Stk Noe paw stk 


We now obtain (4) by substitution of (7) in (6) and by applying 
Lemma 4. In a similar way (5) is obtained by use of Lemma 6. 
Theorem 2 now follows from Theorem 1 by virtue of (3). 


HARVARD UNIVERSITY 


* One may apply, for example, S. Bochner, op. cit., Theorem 12, p. 42. 


— 
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FUNCTIONS OF COPRIME DIVISORS OF INTEGERS 
BY E. T. BELL 


1. Unique Decompositions. If a set U of distinct positive in- 


tegers 1, 1, #2,--~- is such that* 

we call U a coprime set. If to U we adjoin all positive integral 
powers a.>0,--- of integers in U, we 


get the extended set E(U). lf m is in E(U), we call m a U-integer. 


THEOREM 1. Jf n>1 1s representable as a product of powers of 
integers >1 in U, the representation is unique (up to permutations 
of the factors), say 


For, by the definition of U, the u; in (2) are distinct, and by 
(1) a prime p such that p| 7 is such that | u; for precisely one j, 
0<jz<r. We call (2) the U-decomposition of n. 

Obviously there exist U’s such that some 1 >1 are not U-de- 
composable. From the fundamental theorem of arithmetic we 
have the following theorem: 


THEOREM 2. If P=i, pz, --- 1s the set of all positive primes, 
the only U such that every integer n>1 1s U-decomposableis U=P., 


We shall consider also another type of unique decomposition, 
valid for all »>1, which has the distinguishing property of 
U-decomposition as in (2), namely, every n>1 is uniquely a 
product of powers of coprime integers >1. 

If the integer s>0 is divisible by the square of no prime, we 
call s simple. Let S=1, 51, 52, - - - be the set of all distinct simple 
integers; S includes P and is not a coprime set. Without confu- 
sion we may denote by E(S) the set obtained by adjoining to S 


all positive integral powers - - - ,a,>0,a2>0,---, of 
simple integers. 
Let n =p," - - - p,** be the P-decomposition of . If ai, - - - , a, 


are all different, this is by definition also the S-decom position. If 


* In the customary notations, (m, m) is the G.C.D. o. m, n, and m| n signi- 
fies that m divides n arithmetically. 
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are not all different, let a, ---, a; be all the un- 
equal integers among @,---, a,, and let Gig; be all 
those of the a, - - - , a, equal to a;. Writing s;=pa -- - Pig;, we 
have n=s,*!--- si, and this, the S-decomposition of n, is 
unique. 


THEOREM 3. Every integer n>1 is uniquely a product of posi- 
tive integer powers of coprime simple integers >1; the unicity is 
attained when the exponents of the powers are required to be all 
different. 


Note that since E(S) contains E(P), every n>1 has two de- 
compositions, which coincide only if the exponents are all differ- 
ent, into a product of powers of coprime simple numbers, the 
P-decomposition and the S-decomposition as above defined, 
both of which are unique. Thus if 2 = 2?- 3?-53-7-11%, this is the 
P-decomposition, while m = (2-3)?-7- (5-11)? is the S decomposi- 
tion, from the coprime simple integers 2-3, 7, 5-11, with the 
respective exponents 2, 1, 3, all different. 


2. U-divisors, S-divisors. Referring to (2) we define the 

+1) -- - (c,+1) integers 

k, k, : 
(3) U,, 0S k; S ci, 
to be the U-divisors of the n in (2). If m is a U-divisor of n we 
write (m|n)v. Similarly, if 
(4) $= Se, 
is the S-decomposition of s, the S-divisors of s are the (a,+1) 
-- + (a,41) integers 
By an obvious change of notation everything defined next for 
U is defined also for S, and we need state only the definitions 
for U. 

If (m|n)v, there is a unique U-divisor ¢ of m such that mt=n; 
m, t are conjugate U-divisors of n. 

If (d|m)v, (d|n)u, d is a common U-divisor of m, n. If gis a 
common U-divisor of m, n which is such that (d| g)u for every 
common U-divisor d of m, n, then g is unique, and we call 
g=(m, n)y the greatest common U-divisor of m, n. 
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If m, n, t are U-integers such that m=nt, and hence (n| m)v, 
(t| m)v, we call m a U-multiple of n (or of t). The U-integers m, n 
determine a unique U-integer 1={m, n}v, the least common 
U-multiple of m, n, such that if (m|e)uv and (n| e)v, then (i| e)u. 

THEOREM 4. If m, n are U-integers, 


(m, n)u\m, n\y = mn. 


This is proved as in E(P). Let ua, - - - , ua be all the integers 
lay » Ug, Uy, - , usin U occurring in the U-decompositions 
as in (2) of the U-integers m, n, mn~1. Then we may write 

7. he ha 
Ug = Ug Ua, 
t, ts ke ka 


in which some of the h, k may be zero. Writing max(h;, k;) =g;, 
min(h;, k;) =l;, we have 

and hence the theorem. 

If (m, n)y =1, then m, n are called U-coprime. 

In what follows, the particular U-divisors 1m, --- , u, of n as 
in (2) play the part of the distinct primes dividing m in the 
P-decomposition; i, --- , u, will be called the primitive U-di- 
visors of n. By a previous remark, primitive S-divisors are there- 
fore also defined. 


3. Functions of U-divisors. By a change of notation, every- 
thing stated for U holds for P, S, as in §2. The numerous func- 
tions depending on the P-decomposition of integers that occur 
in the theory of numbers,* together with all of their properties 
depending only on the fact that the P-decomposition is unique 
(into a product of powers of coprime integers), go over un- 
changed to the like for U-decompositions by a few obvious 
changes in notation and terminology. We take first the example 
that started the entire theory of such functions for P. 

Let n=p." -- - p,** be the P-decomposition of . Then the 
number ¢() of integers <m and prime to 7 is 


1 1 
(6) = n(1-—)-- 
pi pr 


* See Dickson’s History of the Theory of Numbers, vol. 1, 1919, chapters 5, 
10, 19; also the writer’s Algebraic Arithmetic, 1927. 


Ill 


= 
= 
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The corresponding theorem for U is as follows. Let m=," - - - 
u,“t be the U-decomposition of m. Then the number ¢y(m) of 
integers <m and not divisible by any one of the primitive U- 
divisors %,--- ,u,of mis 


(7) = ~) (1 - ~), gu(1) = 1. 
My 


The proof of (7) is precisely similar to that of (6) by means of 
the principle of cross-classification,* with the remark that 
(a, 6) =1, b| together imply ab] k. From the explicit form 
of ¢() in (6), a common algebraic proof gives Gauss’ result 


(8) > o(d) = n, 
Hence (7) implies 
(9) ou(t) =m, (t| 


Generally, to pass from P to U, PU, we have 
U, 

“prime” — “primitive”, 

“divisor” — “U-divisor”, 

(m,n) —(m, n)v. 


(10) 


If f(x) is single-valued and finite for integer values >0 of x, 
f(x) is (as usual) called a numerical function of x. The unit nu- 
merical function n(x) is defined by =1, n(x) =0, x¥1. A 
generalization of Dedekind’s inversion formula, proved in a pre- 
vious paper, is of great use in the algebra of numerical func- 
tions. If, and only if, f(1) +0 there exists a unique numerical 
function f’(x), such that 


(11) > f(d)f’() n(n), 1, 2; | 


the sum referring to all pairs d, 6 of conjugate P-divisors of n. 
Passing from P to U by means of (10), we get the theorem cor- 
responding to (11) on making the single change U for P in the 


* First stated as a general principle in arithmetic, apparently, by da Silva 
in 1854, Memorias da Academia Real das Sciencias de Lisboa, N.S.I., pp. 8-9; 
see Dickson, loc. cit., p. 119. A special form of the principle was noted in 1857 
by H. J. S. Smith; see his Collected Mathematical Papers, vol. 1, p. 36. 

t Téhoku Mathematical Journal, vol. 17 (1920), pp. 221-231. Simplified 
proof, ibid., vol. 43 (1937), pp. 77-78, 


= 
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foregoing statement. The proof may be given precisely as in 
the references cited. The generalization mentioned is as follows. 
If f(x), g(x), h(x) are numerical functions such that 


> fags) = h(n), n=1,2,---, g(1) 
> g(d)g'(6) = n(n), 


ll 
N 


then 
(12) f(n) = 


all sums referring to all pairs d, 6 of conjugate P-divisors of n. 
If g(n) =u(n), =1 for n=1, 2,---, g’ is Mobius’ p, and (12) 
becomes Dedekind’s inversion. To pass from (12) to its U-corre- 
spondent it suffices to replace P by U as for (11). The U-corre- 
spondent of Dedekind’s inversion is obtained by replacing 
“conjugate P-divisors” by “conjugate U-divisors,” and yp by yu, 
where uy (n) is zero if is divisible by the square of any primitive 
U-integer, and otherwise is +1 or —1 according as 7 is the 
product of an even or an odd number of primitive U-divisors 
of n; by convention py(1)=1. A similar convention holds for 
any fu() which is not otherwise defined when n=1, namely, 
fu(A) =1. 

In previous papers* the algebra of numerical functions based 
on P-decomposition was constructed from (11), (12). It follows 
that there is a simply isomorphic aigebra for any U-decomposi- 
tion. In the P-algebra it was noted that the theorems hold for 
any set in which there is a unique decomposition. Hence the like 
is true for U-decompositions in such a set, for example the ideals 
of an algebraic number field. 


CALIFORNIA INSTITUTE OF TECHNOLOGY 


* Some are listed in my paper, Journal of the Indian Mathematical Society, 
vol. 17 (1927-28), pp. 249-260. 


1937-] NEO-SYLVESTER CONTRACTIONS 823 


NEO-SYLVESTER CONTRACTIONS AND THE 
SOLUTION OF SYSTEMS OF LINEAR 
EQUATIONS* 


BY F. L. WREN 


1. Introduction. In 1851 Sylvester{ stated, without proof, a 
theorem on contraction of determinants which had previously 
been stated by Hermite (1849)f{ for a special case. A proof of 
Sylvester’s theorem was given by Studnitka in 1879.§ A re- 
statement of this theorem in a slightly different form permits 
significant applications to the determination of the ranks of 
matrices, the solution of systems of linear equations, and the 
calculation of partial and multiple coefficients of correlation. 


2. Restatement of Sylvester’s Theorem. The restatement of the 
theorem is as follows: 


THEOREM 1. Jn the nth order determinant, 


D= |ai;|, (i,j =1,2,---,), 
let 
kt, ° °° 
= 


be a non-vanishing minor of order r, and associate with each ele- 
ment a;; (t%#k, 741) an (r+1)-rowed minor of D defined by the 
identity 


* Presented to the Society under another title, December 29, 1936. 

1 J. J. Sylvester, On the relation between the minor determinants of linearly 
equivalent quadratic functions, Philosophical Magazine. (4), vol. 1 (1851), pp. 
295-305, 415. 

¢ C. Hermite, Sur une question relative a la théorie des nombres, Journal de 
Mathématiques, vol. 14 (1849), pp. 21-30. 

§ F. J. Studnitka, Ueber eine neue Determinantentransformation, Sitzungs- 
berichte Geschichte der Wissenschaften, vol. 9 (1879), pp. 487-494. 


A 
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aij ail, 
Dij =| * * 
kt, * * k,l, 


in which M,, ts the first minor corresponding to a;;. We have, then, 
the relation 


n—r—l 


(— 1)* My D=A, 
+h, 


where A is a determinant of order (n—r) defined by the identity 


Dii, 
A= GARD, 
in which the (t, %t2,---, tn-r) represent the rows and 
(ji, je, » Jn—r) the columns of the minor of D complementary 


to 


CoROLLARY. Attach to the determinant D a row whose elements 
are the sums of the elements in each respective column of D, and a 
column whose elements are the sums of the elements in each respec- 
tive row of D. If the elements of the added row and column be 
treated in the same manner as the elements a;;, (i~k, 71), were 
treated in the theorem, then the values obtained will be the sums of 
the elements D;; in the corresponding column or row of A. 


The determinant D with the attached row and column is 


* For a proof of this theorem see G. Kowalewski, Einfiihrung in die De- 
lterminantentheorie, pp. 90-93. Kowalewski points out that the theorem holds 
true even for M;,=0. Attention is also called to the fact that when the minor 
Mi: is a single element the expansion is the same as that given by Chié (1853) 
in his Mémoire sur les foncticns connues sous le nom de résultantes ou de détermi- 
nants. 
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5; 


n 
a1 12 aij 
j=1 


n 


(1) j=l 


n 


n n n 
i=1 i=1 i=1 


If the elements in the S,-column and the S,-row are treated as 
the elements a;;, then the resulting A, with its corresponding 
S,-column and S,-row, is 


Diu Diz Din | 
Day Dee Den 1) 


where the D;; are defined for i+k and j#/, and 


n 

> aij il, 

j=1 

n 

Gk yj Gk ak ak 
D;; j=1 


n 
j=1 


= 
= 
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By subtracting from the first column of D,; the sum of the re- 


maining columns, we have 


= 
i 


where the summation on j is extended over the numbers 


A similar argument gives 


, n, excluding lo, - - - 


D.; Di;, 


where the summation on 7 is over the numbers 1, 2, - - - 
- , k,. (It might be added that, for all practical 


cluding ki, ke, 


, N, ex- 


purposes, either the S, row or the S, column is sufficient in itself 
for a satisfactory relative check.) 


3. The Rank of a Matrix. If a matrix is not of zero rank there 


is at least one element a;;+0. 


If, in the mXn matrix 
fair 


d21 


Ami 


G@i2 °° Gin 
d22 den 
Cnn) 


of rank +0, an element a;;+0, called the pivotal element, be 
chosen and 2 new matrix be formed whose elements are the sec- 


ond order determinants 


then a new matrix 
Du 
M’ = | Ps 


Dm 


ai 
Diz Dis 
Deo Den 
| 
D2 


will be formed. This matrix will have (m—1) rows and (n—1) 
columns, since there are no elements in the kth row and the /th 
column. Furthermore, the elements D;; are all of the possible 


non-vanishing 2-rowed minors of M of which a,, is a first minor. 


= 
= 
a; 
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Hence if all the D;;=0, the matrix M is of rank 1. If there is an 
element D,,~0, proceed as in the previous case, using D,, as the 
pivotal element to construct the (m—2) X(n—2) matrix 


(Arn Ars Ain 


= | Aor Ase |, 
| 


(Ami Ane Ane ) 


where the elements are as follows: 


| Dypq Dys Qpq | 
Get | eh 
Qkq Akl Ges Akl 
Hence, by Theorem 1, 
Ips pt | 


k, ‘) 
Ang = Ors Ort |, 
q#l,s 
Thus it is seen that the determinants 4 ,, are, except for the non- 
vanishing factor a,;, the possible non-vanishing third-order mi- 
nors of the original matrix M of which the non-vanishing 
second-order minor D,, is a first minor. Hence, if all A,,=0, the 
rank of M is 2. 

If there is a A,,#0 the process may be repeated, this time 
using the elements in M’’. This application, by use of Theorem 
1, will give a matrix whose elements will be 

a41-D,,X (all possible non-vanishing third order minors of M 

of which D,, is a first minor). 

Repeated applications of the above contraction process, 
which will be called a neo-Sylvester contraction, give the follow- 
ing result: 


THEOREM 2. A matrix is of rank r if r neo-Sylvester contractions 
reduce it to a zero matrix. 


CoROLLARY. If, in an m Xn matrix (mn), the result of m 
neo-Sylvester contractions is a non-zero matrix, the rank of the 
original matrix is m. 
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The check established in the Corollary to Theorem 1 also ap- 
plies in the treatment of matrices. 


4. Systems of Simultaneous Linear Equations. The augmented 
matrix of a system of m linear equations in unknowns is the 
m X (n-+1) matrix 


| a21 d22 dan ke 
Amn 


Since the coefficients a;; appear simultaneously in the deter- 
minants of both the coefficient and augmented matrices, it is 
obvious that the ranks of these two matrices can be obtained 
simultaneously by selecting for each contraction a pivotal ele- 
ment from the coefficient matrix. 


THEOREM 3. A system of m linear equations in n unknowns is 
inconsistent if and only if a series of neo-Sylvester contractions, 
with pivotal elements chosen each time from the coefficient matrix, 
will reduce the augmented matrix to a matrix with zeros everywhere 
except in the column of known terms k;. 


Each such contraction of the matrix of the system of linear 
equations is €quivalent to eliminating the unknown that corre- 
sponds to the column of coefficients from which each successive 
pivotal element is chosen. Consequently, when the system of 
equations is a consistent system of rank r, the result of (r—1) 
neo-Sylvester contractions is the non-zero matrix of a system of 
(m—r+1) equations in (7 —7r+1) unknowns in which the corre- 
sponding coefficients are proportional. Any one of these equa- 
tions may be solved for one of the (7 —r+1) unknowns in terms 
of the remaining (n—r). This value may be substituted in one 
of the equations obtained by (r—2) contractions. Continued 
substitution in the next previous system of equations results 
in the final determination of r of the unknowns in terms of the 
remaining (n—r) unknowns. These values will satisfy the re- 
maining (m—r) equations of the resultant (m—r+1) linearly 
dependent system obtained by the (r —1) neo-Sylvester contrac- 
tions. 


A 
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THEOREM 4. If (r—1) neo-Sylvester contractions reduce the 
augmented matrix of a system of m equations in n unknowns to 
a non-zero matrix of dimensions (m—r+1)X(n—r+2) in which 
the elements of any one row are proportional to the corresponding 
elements of every other row, then r of the unknowns can be deter- 
mined uniquely in terms of the remaining n—r unknowns. These 
values will also satisfy the remaining m—r equations. 


If r=m =n, we obtain the following result: 


THEOREM 5. If (n—1) neo-Sylvester contractions reduce the 
augmented matrix of a system of n equations in n unknowns to a 
1X2 non-zero matrix, then the equations have a unique solution.* 


Each contraction is equivalent to the elimination of variable 
x; if the pivotal element is chosen in the jth column. Hence, if 
the pivotal elements for the (n—1) contractions be chosen from 
columns other than the column of coefficients of x, and the col- 
umn of constant terms, the value for x, is obtained by dividing 
the final element in the column of constant terms by the final 
element in the column of coefficients of x,. This is obvious, since 
the final element in the x, column is, except for a sign and multi- 
plicative factor, the expansion of the determinant of coeffi- 
cients, and the element in the column of known terms ; is, 
except for the same sign and multiplicative factor, the expansion 
of the determinant obtained from the determinant of the coeffi- 
cients by replacing the elements of the rth column by the col- 
umn of k’s. 


Corotuary. In the nX(n+1) matrix 


(dir Gin | 
G22°°* Gen 
M = 
Qn1 Gn2°** Ann 


let (n—1) neo-Sylvester contractions be made by choosing the piv- 
otal element for each contraction from any column other than the 
jth or kth columns. The quotient of the final element in the kth 


* A numerical example of the truth of this theorem and the use of the check 
suggested in connection with Theorem 1 is found in H. G. Deming’s article, 
A systematic method for the solution of simultaneous linear equations, American 
Mathematical Monthly, vol. 35 (1928), pp. 360-363. 


— 
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column by the final element in the jth column is, except for sign, 
the quotient of the two determinants of order n, 


G21 422 den a2; den 


in which the elements are all alike except those of the jth and kth 
columns. The sign factor is (—1)!*-‘!-! where |k—j| —1 ts the 
number of inversions necessary to place the kth column in the divi- 
dend in the same relative position as the jth column in the divisor. 


The truth of this corollary is readily seen. The contractions 
would produce in the jth column of M the expansion of the de- 
terminant which is the divisor in the Corollary except for a mul- 
tiplicative factor, which would be the product of the (m —1) piv- 
otal elements chosen for the successive contractions, and for a 
possible sign factor. The same contractions would produce in 
the kth column of M, except for the same sign and multiplica- 
tive factor, the expansion of the determinant which is the divi- 
dend in the Corollary multiplied by the additional sign factor 
(—1)!4-4I-!, which is the number of inversions necessary to 
place the kth column of the dividend in the same relative posi- 
tion to the remaining columns as the jth column of the divisor 
bears to these same columns. 


5. Partial Correlation. In the n-variable correlation problem 
the regression equations are 


(i, 7=1, 2,---,m except that neither has the value 9p), 


where r;; are ordinary product-moment coefficients of correla- 
tion between variable i and j (hence 7,;=r;; and r;;=1) and 
Bpj 12---(p)+--G)+--n are the regression coefficients that minimize 
the error in estimating variable p from a linear function of the 
remaining (7—1) variables.* These are (n—1) linear equations 
in (n—1) unknowns and hence the use of neo-Sylvester con- 


* Truman L. Kelly, Partial and multiple correlation, in H. L. Rietz, Hand- 
book of Mathematical Statistics, 1924, pp. 139-141. 


| 
| 
| 
| 
| 
= 
| 
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tractions affords a systematic technique for determining the val- 
ues of the 6’s; furthermore the check used in solving linear equa- 
tions may be used here. 

The process may be further systematized by making use of 
the mth order determinant of the product-moment coefficients 


A= 


where A,, will represent the minor of order (nm — 1) obtained from 
A by deleting the pth row and qth column. 

It is to be noted that the determinant of the coefficients of 
the 6’s in the regression equations (2) is the determinant A,, and 
that the known terms are elements from the pth column of A. 
Perform (n—2) neo-Sylvester contractions on A in which the 
pivotal element is chosen each time as any element in the prin- 
cipal diagonal except in the pth or gth row and column. The re- 
sult of these (7 —2) contractions will be 


p-column g-column 
(3) p-row | KAg, (— 
g-row | (— 


where K is a common sign and multiplicative factor. Since A is 
symmetric and the pivotal element for each contraction is 
chosen from the principal diagonal of A, it follows that A,, 
=A,,; furthermore the sign factor may be written (6+q—1), 
since, for all integral values of p and g, (b-+q—1) is even or odd 
according as | p—q| —1 is even or odd. Hence the determinant 
(3) may be written in the simpler form 


p-column g-column 
(4) p-row | (— 1)***'KA,, 
g-row | (— KA,,| 


Therefore, by the Corollary to Theorem 5, 

(= (= 1)? 

KAgq 


(5) 


| 
| 

\ 


832 F. L. WREN (December, 


The values of the remaining B);.12...;p)..-)---n Can be de- 
termined as in §4. Furthermore, if the first pivotal element is 
chosen in the kth row and column the following 6’s may be ob- 
tained as described : 

B;;..: Divide the element in the ith row and jth column by 

the element in the jth row and column. 

B;:-.: Use the process symmetrical to the one above, that is, 
divide the element in the jth row and ith column by the 
element in the ith row and column. 

If the second pivotal element is chosen from the /th row and 
column, §;;.4: and 8;;... may be obtained as above. Each subse- 
quent contraction will produce the conjugate regression coeffi- 
cients of i on j and j on 7 with the additional variables controlled 
corresponding to the row and column in which each subsequent 
pivotal element is chosen. 

The partial coefficient of correlation fpq.12...¢p)---(q)---n has 
the value (—1)?**"1A,,/(AppAqq)’?.* Hence the value of 
1 pq-12---(p)-+-(q)---n May be obtained directly from determinant 
(4). Furthermore, 7;;., may be obtained after the first contrac- 
tion using the element in the kth row and column as the pivotal 
element. It is obtained, as above, from the two rowed minor of 
the contracted determinant obtained by deleting all rows and 
columns except the ith and jth. Similarly, after each subse- 
quent contraction, values may be obtained for 


Vij-kly Vij-klmy Tij-klmny * * 


Furthermore, each contraction can be checked by attaching to A 
a column whose elements are the sums of the elements of each 
respective row of A and applying the same contraction to these 
elements that is applied to the elements of A. 


6. Multiple Correlation. The multiple alienation coefficient is 
given by 


2 A ft 
App 


If we define the determinant A (p) by 


* Truman L. Kelley, Statistical Methods, p. 298. 
t Ibid., p. 301. 
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fis 1 
Tip 3p 1 pn 


where the elements of the pth column of A have been replaced 
by zeros except for the element 1 in the pth row and column, 
it is evident that 


A(p) = App. 
Formula (6) may now be written 
(7) 
A(p) 


If the determinant of the numerator be augmented by adding 
the pth column of the denominator, the »X(n+1) matrix 


(1 Tis Tip Tin | 

(8) | Sip Tap 1 1 | 


is obtained. If this added column be called the p’-columr of the 
matrix, and if the pivotal elements for the neo-Sylvester con- 
tractions be chosen each time from columns other than the p- 
and p’-columns, then by the Corollary to Theorem 5, after 
(n—1) such contractions the value of k?.12...;»).--» can be ob- 
tained by dividing the element in the p-column by the element 
in the p’-column. Furthermore, the result of each contraction 
can be checked by the method of the Corollary to Theorem 1. 
The multiple correlation coefficient is given by the relation 


2 2 
(9) Tp-12---(p)---n = 1 — 


z 
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Substituting in (9) from (7) and simplifying, we have 
_ AG) 


Tp-12---(p)---n = 


A(?) 
1 T12 «713° Tis Tin 
(10) 
Tin Ten 1 


A(p) 


Thus, by augmenting the numerator of (10) as in the case of 
(7), the pth column of the denominator becomes the p’th column 
of the matrix 


(1 13 °** te 

rio 123 — Tep fs, 0 
(11) 


Hence the same process and check which was applied to matrix 
(8) to find R?.12.. .;»)---n, When applied to matrix (11), will give 
1; .12---(p)---n- Furthermore, if the rows and columns containing 
the pivotal elements are chosen in the sequence 1, 2, 3,---,m 
(omitting p) and if, after each contraction, the element in the 
pth row and column be divided by the element in the pth row 


and ’th column, then the values k?.;, k?-12, k?-123, , and 
72 4, 12.12, 72.123, -- - Will be obtained from matrices (8) and (11) 
respectively. 
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NOTE ON AN ELEMENTARY GEOMETRIC 
EXISTENCE THEOREM 


BY P. W. KETCHUM 


1. Introduction. By the []-moment of a set of points 
C1, C2,°-* , Cm, (m>1), with respect to a particular c, we mean 
the product of the m—1 distances from c, to the other c’s. The 
object of the present note is to prove the following proposition: 


THEOREM 1. For any given set of distinct points ay, d2,--- , Gn, 
(n>1), in the plane, there will exist a set of points bi, be, --- , On, 
distinct from the a’s, such that the | |-moment of the a’s and b’s with 
respect to a; is the same for all i. 


This theorem is obviously a special case of the following more 
general theorem: 


THEOREM 2. Let a, d2,--- ,@n, (n>1), be a set of distinct 
points in a plane and let d;,; be the distance from a; to a;, (t#j).* 
Let r; be any given ray (half line) issuing from a;. Then for every 
sufficiently small positive number 6, say for 5S5o, there exists a 
set of points by, be, --- , b, satisfying the two conditions (a) the 
point b; is on the ray r;, and (b) if D;,; 1s the distance from a; to b;, 
then 


for all values of 1. 


2. An Outline of the Procedure. An algebraic proof of the exist- 
ence of the desired solutions of (1) seems to be difficult, so we 
prefer to use the method of successive approximations. The 
formal details of our procedure are a little tedious, but the un- 
derlying idea is relatively simple, namely: Since 6 is small, the 
point b; will be much closer to a, then to any of the other a’s or 
b’s. Thus, if b; is shifted slightly, the quantity D;,; will change 
relatively much more than any of the other distances. This has 
the effect of making the product, [], in (1) approximately a 


* Throughout this note, 7 and j will range over the values 1, 2, 3,---, n, 
and k over the values 1, 2, 3, - - - 


= 
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function of one variable D;,; only. As a first approximation we 
therefore assume that D; ;=d;,; for 1~j, and calculate D;,; so 
that the product [| will equal 5. This will determine new dis- 
tances D;_;, i*j, slightly different from d;,;, which one may use 
as a second approximation for the calculation of D;,;; and so on. 
It is intuitively evident that the process will be a convergent 
one. We now proceed to the details. 


3. Proof of Theorem 2. We are going to prove the theorem 
for the particular value 


59 = where D = max (d;,;), d = min (d;,;). 


Put D,9 =d; ; for ixj, and let D,) be a set of numbers satis- 
fying the relations 


(2) 
(i) (k~-1) 


II (d;,;Di,; ) 


Here [[,“® is a symbol indicating a product to be taken over all 
values of 7 with the exception of the value 7=17, which is to be 
omitted. Denote by b;*) the point on 7; at distance D,* from 
a;. Let DP be the distance from a; to b;*). The quantities 
D{? and b{* are thus uniquely defined by recurrence. 

We first show that the points 5;*) have for each 7 a unique 
limiting point as ko. Set 


(1) (1) (1) (1) (0) 


nii = = for 
and 


(k+1) (k+1) (k) 
ng = — Di,;- 


It will then be sufficient to prove that 


(k+1) —1_2—5n—3k n 


(3) | <n 2 


since the right-hand member is the kth term of a convergent 
geometric series. We proceed to prove (3) by induction. Put 


(k) (k) 
ne = max | 9:,;| = max | 7;,;|, 


Do = Dy-1 + do = d, = 


= 
= 
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Then 
(k) (k) 


Di,; Dz, D;,; 2 dx. 
Also 


(k) (k—-1) 


(i) 


Daal 


(k4+ +1) 
Nii 


Calling the numerator in this expression 6K*+» and the de- 
nominator 7 we have 


_(k+1) | 


| lot 1) "| 


IIA 


n—2 n—3 n—1 
(n = + ; 
where the C’s are binomial coefficients; and 


(k+1) n—1 n—1 n—1 


T; 2d dy_1dy 


provided d; is positive for all k. Now from (2) and our hypothe- 
sis on 6, 


3 
ll 
t 


Hence 


= 
IIA 


An easy computation shows that 


(2) —1_—3n .n—1 


and 
(2) 
nil = <n 2 D 


Thus (3) is seen to be satisfied for the first value of k, namely 
k=1. 


| 
Now 
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Next we assume that (3) is satisfied for k=1, 2,---,u—1. 
Then 


An easy computation shows that 


Now 

D, =D+m+2+--- +m 

1 1 
| =2D, 
from our assumption that (3) holds for k=1, 2,---, w—1. 
Hence 
< 

Also 

d,=d—m—m2- — 

1 1 
8 16 
Hence 
re 
and 
ett) 
T;@t) 


so that (3) is true when k=y. The induction is thus complete 
and (3) is proved, and the points };*) have a unique limit point. 

Now take the points 6; of Theorem II as the limiting points 
of 5. Let 


(i) 


(k) (k) 
6; = II 4:.; I] D;j,;. 
7 
Since : 
D;,; = lim 


kx 
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we have 
(i) 


ij i 


Now 


But the quantity on the right approaches zero, so that 5{*) -6 
as k->«. We thus have (1), and the theorem is proved. 


INSTITUTE FOR ADVANCED STUDY 


AN INVOLUTORIAL LINE TRANSFORMATION IN S, 
BY C. R. WYLIE, JR. 


1. Introduction. It is a well known fact that all planes which 
meet four general lines of S,;* are met by a fifth line. The re- 
markable configuration determined by five such “associated 
lines” is discussed in a number of places in the literature.f In 
the present paper an involutorial line transformation suggested 
by the figure of five associated lines is discussed, both as a line 
involution in S,, and as a point involution on a certain V¢ in So. 
In §§2-6 the involution is treated at some length by purely syn- 
thetic methods. The final section (§7) contains a brief analytic 
treatment, including the equations of the involution, and the 
equations of the invariant and singular elements. The involu- 


* We shall use the conventional symbol S,, to indicate a linear space of 
dimension m. A variety of order r and of dimension m we shall designate by the 
symbol! 

{t Welchman, W. G., Plane congruences of the second order in space of four 
dimensions and fifth incidence theorems, Proceedings of the Cambridge Phil- 
osophical Society, vol. 28 (1931-1932), pp. 275-284. 

Baker, H. F., On a proof of the theorem of a double six of lines by projection 
from four dimensions, Proceedings of the Cambridge Philosophical Society, 
vol. 20 (1920-1921), pp. 133-144. 

Baker, H. F., Principles of Geometry, Cambridge University Press, 1925, 
yol. IV, Chapter V, ; 


(i) 
(k— 
II 
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tion appears to be an important one in that most of the loci con- 
nected with it are of considerable interest in the geometry of 


hyperspace. 


2. Representation of Five Associated Lines in Sy. Stephanos 
has shown* that five associated lines have their ten Grassmann 
coordinates linearly related. If we regard the ten Grassmann 
coordinates of the lines of S,; as point coordinates in S»,f and re- 
call that the five quadratic identities existing among these co- 
ordinates define a V,° in Syf whose points are in 1:1 correspond- 
ence with the lines of S,, this result may be interpreted as fol- 
lows: The images in Sy of five associated lines of S, are the five 
points in which a general S; meets V5. In fact, from this as a 
definition of five associated lines the incidence property men- 
tioned above follows at once. For consider four general lines, 
L,, Le, Lz, Ls, of Ss, and any plane, P, meeting each of these 
lines. On V,° we have the four independent points Js, Ls, 
which are the images of the four lines Z;, Lz, L3, L4.§ Now the 
Ss which intersects V,° in the V;> representing the special com- 
plex of lines which meet the plane P contains |, h, Js, 1. It 
therefore contains /;, the fifth point in which the S; determined 
by hy, le, ls, 1, intersects V,°. Hence /; is also the image of a line. 
L;, which meets the plane P. 


3. Definition of the Involution. Let there be given three gen- 
eral lines, Z;, Le, L3, of S,; and let the transform of any line L 
be the line LZ’ which forms with LZ, Le, Z3, and L a set of five 
associated lines. On V,° this becomes a point involution defined 
as follows: Let there be given three points, /), /2, /;, on Ve> such 
that the plane, 7, which they determine meets V,° only in these 
three points. Then the transform of any point, /, is the fifth 
point, /’, in which the S; determined by h, /2, ;, and / meets V,°. 


* Stephanos, C., Sur une configuration remarquable de cercles dans l'espace, 
Comptes Rendus, vol. 93 (1881), p. 578. 

+ Todd, J. A., The locus representing the lines of four dimensional space 
and its application to linear complexes in four dimensions, Proceedings of the 
London Mathematical Society, (2), vol. 30 (1929-1930), pp. 513-550. 

t See §7 for the analytic definition of the Grassmann coordinates, and of 
the five quadratic identities. 

§ Throughout this paper we shall represent configurations in S, by capital 
letters, and their images in Sy by small letters, 
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To determine the order of the involution consider a general 
pencil of lines, R, of S,;. The image of this pencil on V,° is a line, r. 
The «' spaces S; determined by z and the respective points of r 
all lie in the S, determined by z and r. This S, meets V,° in a 
curve of order five. The line 7 is a point of this curve of inter- 
section. The other component is a rational normal quartic 
curve, passing through the three points /;, 2, /; in r. This quartic 
curve is the transform of r in the involution on V,°, hence the 
involution in S,; is of order four, transforming a pencil of lines 
into a ruled quartic surface, three of whose generators are the 
fundamental lines L;, Le, L3. 


4. Invariant Elements. There are ~‘ varieties V? on V,°, these 
being the representations in Sy of the lines of the 4 spaces S; 
in S,.* Each of these V?’s lies in an S;. One and only one of these 
S;’s passes through a general point of Ss, while »? of them pass 
through each point of V,*. In particular the ~? spaces S; corre- 
sponding to the points of 7 cut V,> in V?’s each of which contains 
the point, ¢, which is the representation in Sy of the common 
transversal, 7, of the lines Z;, Le, Z3, in Sy. These V?’s are 
thus representations of the lines of the ~* spaces S; through T 
ye 

Consider now « point, /, and its transform, /’, in the involution 
on V,°, and consider further the S; corresponding to the point, q, 
in which the line joining / and /’ intersects 7. Any other point, m, 
on the V? cut from V,° by this S; must have for its transform, 
m’, another point of the same V?, namely the second point in 
which the line joining g and m intersects the V?. Thus the in- 
volution on V,° leaves invariant as a whole, though not point by 
point, the ~? varieties V? corresponding to the points of 7. This 
implies that in S, the transform of any line L lies in the S; de- 
termined by Z and the transversal, 7, of the three lines Zi, 
Lo, Ls. 

The invariant points in the involution on V,° are the points 
of contact of the tangents to V,° which meet z. A triple infinity 
of these tangents pass through each point, g, of 7. They are in 
fact the »* tangents to the V? corresponding to the point gq; 
and the locus of their points of contact is the V? cut from the 
V2 by the polar S, of g as to the V? (the polar S, being deter- 


* Todd, loc. cit. 


> 
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mined in the S; containing g and the corresponding V?). Thus 
in S, the invariant lines of the involution form a linear complex 
in each S; through T. 

On V,° there are ~?* planes with the property that they con- 
tain three non-concurrent lines which determine with the respec- 
tive points /,, /2, 1; planes lying entirely on V.°. (These planes are 
the representations in Sy of the fields of lines lying in planes 
which meet Z, Le, L3.) Each of the three lines in such a plane, ¢, 
is singular in the involution, its image being the plane which 
joins it to the corresponding point, /;, plus two lines joining /; to 
the other two points, /; and /,. The image of ¢ itself is the resid- 
ual V+ cut from V,> by the S; determined by z and oc. This V# 
is composite, consisting of the three planes which are the images 
of the three singular lines in ¢, and a fourth plane, o’. Now oa’ 
must meet in a line each of the other planes composing V2, 
since it lies on V,° and meets in a line the S,; determined by 7+ 
and each of these planes. But this requires either that o’ coin- 
cide with o, or that o and the planes composing V+ lie in an S, 
instead of in an S;. The latter is impossible, hence o and a’ co- 
incide, and the non-singular points of o are invariant. Con- 
versely, every invariant point on V,° lies in a plane o, namely the 
plane common to the three V?’s determined by the invariant 
point and h, respectively. 

Thus in S, the invariant lines are the ~5 lines which lie in 
planes meeting Z;, Le, L3. Hence in a general S; of S, the in- 
variant lines form a tetrahedral complex, whose fundamental 
tetrahedron is determined by the four points in which the 5; in- 
tersects L,, Le. L3, and 7.* In an S; through 7 this quadratic 
complex breaks up into a general linear complex, as we have 
already noted, and the special linear complex of lines meeting 7, 
the latter complex being singular and not properly invariant. 


5. Singular Elements. The singular elements in the involution 
on V,°, aside from the three fundamental points h, /:, 13, are 
those points which determine with 7 an S; meeting V,° in a 
curve instead of in five points. They are of three general types: 
Those points which lie on lines of V.° passing through /;, those 


* The general tetrahedral complex can be formed by sectioning with an S; 
the system of all planes meeting three general lines of S;. Cf. Baker, Principles 
of Geometry, Vol. IV, p. 32. 
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points which lie on conics of V,> passing through /; and /;, and 
those points which lie on cubics of V¢° passing through /;,1;, 1x. 
The last class of points is also the class of points lying on lines 
of V5 which pass through ¢. This is evident when we consider 
the S, determined by ¢ and an S; containing a cubic curve 
through /,, 2, ;. Such an S, contains not only the cubic curve 
but also the lines joining ¢ to h, h, Js. Since these lines lie on V,°, 
the S, thus meets V,° in a curve of order greater than five, and 
hence has a surface in common with V,°. This surface must be a 
cubic cone with vertex at ¢, and the given cubic as directrix 
curve. Hence all points on any cubic through h, h, J; lie on lines 
through 

The three classes of singular points correspond respectively 
to the following classes of singular lines in S,: Lines meeting L;, 
lines lying in the S; determined by L;and L;, and lines meeting 
T. The images of lines belonging to one or more of these classes 
may be described as follows: (The symbol ~ is used to mean 
“is transformed into.” ) 


1. A general line, L, meeting 7~ the cubic regulus containing 
I, L2, L3, and L as generators, and having T as directrix. 

2. A general line, Z, in the S; determined by L;and L; ~ the 
quadratic regulus determined by L;, L;, and L. 

3. A general line, L, meeting ZL; ~ the pencil determined by 
and L. 

4. A general line, Z, through the intersection of T and L; ~ 
the pencil determined by T and L,, together with the quadratic 
regulus determined by L;, L;, and the line of the pencil which 
lies in the S; determined by L; and Ly. 

4.1. A general line, ZL, meeting T and lying in the 5S; deter- 
mined by L;and L; ~ the quadratic regulus determined by L;, 
L;, and L, together with the pencil determined by L; and the 
generator of the regulus which passes through the intersection 
of 7 and L,. 

5. A line, LZ, in the S; determined by L; and L; passing 
through the intersection of T and L; ~ the pencil determined by 
L, and L, the pencil determined by L; and the line of the first 
pencil which meets L;, and the pencil determined by L; and 
the line of the second pencil which meets Ly. 

5.1. A general line. L, in the plane of T and L; ~ the pencil 
determined by L; and L, the pencil determined by L; and the 
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line of the first pencil which meets L;, and the pencil determined 
by L, and the line of the first pencil which meets L;. 

6. A general line, L, meeting L; and L; ~ the pencil deter- 
mined by L; and L, together with the pencil determined by L; 
and L. 

6.1. A general line, L, meeting L; and lying in the S; deter- 
mined by L; and L; ~ the pencil determined by L; and ZL, to- 
gether with the pencil determined by ZL; and the line of the first 
pencil which meets L;. 

7. The transversal, 7, ~ the three pencils determined respec- 
tively by JT and L;, T and Lj, and T and Ly. 

8. The fundamental line, L;, ~ the special linear complex 
consisting of all lines which meet the plane determined by T 
and L,. 


None of these but the last needs special comment; we verify it 
in this fashion: Any point which is an image of /; in the involu- 
tion on V,° must be such a point that the S; which it determines 
with x contains one of the tangents to V,° at /;. There are ~° 
tangents which can be drawn to V,° at the point /;. All these 
points lie in the Ss determined by 7 and the tangent S, to V,° 
at /;, and hence correspond to lines of a special linear complex 
in S,. The singular plane of this complex must contain L;, and 
must be met by L;, and L,. It must therefore be the plane de- 
termined by 7 and L,. 


6. Images of Linear Systems of Lines. As we have already 
noted, every line, Z, which meets either Zi, Z2, or L3, is trans- 
formed into a pencil containing L. Although such lines are singu- 
lar they may thus be regarded as invariant also. Hence Zi, Le, 
and ZL; must lie upon all the quadratic hypercones formed by 
the invariant lines which pass through the respective points of 
S;. Moreover L;, Le, and L; lie in planes of the same family on 
each of these hypercones. 

Consider now the ~# lines through a general point of S,. On 
V° these are represented by the points of an S;. The transform 
of this S; is the residual V34 in which the S, determined by 7+ 
and the 5S; intersects V.>. Now the Ss whose intersection with V¢° 
represents the special complex of lines meeting any plane of the 
second family (the family not containing L;, L2, and L3) on the 
hypercone of invariant lines through the given point in S,, con- 


1937-1 LINE TRANSFORMATION IN S, 845 


tains ],, , ; and the S; which represents the totality of lines 
through the given point. The Ss therefore contains the entire S¢ 
in which the image V3 lies; hence all points of V/ represent 
lines of S, which meet all planes of the second family on the 
hypercone of invariant lines. Thus V;‘ must represent the total- 
ity of lines lying in the planes of the first family on the hyper- 
cone, and this system is then the transform of the ~# lines 
through a general point of S;. 

Since the double infinity of lines through a general point in a 
general S; in S, is included in the totality of all lines through the 
point, the transform of such a system must be ~? lines lying in 
planes of the first family on the hypercone of invariant lines 
through the point. These lines are in fact the lines which lie in 
planes of the first family and meet the plane determined by the 
points where the 5S; intersects L;, L2, and L3. For on V,° the 
image of the given system of lines is a plane, and its transform 
is the residual V3 in which the S; determined by z and this 
plane meets V,°. Now the Ss which meets V,° in the V;° repre- 
senting the special complex of lines meeting the plane, P, de- 
termined in S, by the points where the given S; intersects Li, Lo, 
and L3, contains /;, 2, /;, and the plane which represents the 
given system of lines. It therefore contains the S; in which the 
image V+! lies. Hence all points of this Vs! represent lines meet- 
ing the plane P. Since the family of lines represented by this 
V+ contains «! penci's of lines (one in each plane of the first 
family on the hypercone), the V2 must be the rational normal 
ruled V2! of S;, and not the Veronese surface. 

We have already noted that the transform of a general pencil 
of lines of S, is a rational ruled quartic surface containing Ly, Le, 
L;, and two lines of the given pencil, namely the two invariant 
lines. These two lines are the only generators of the surface 
which can intersect. This surface, the projection of the rational 
normal ruled quartic surface of S;, must also lie on the hyper- 
cone of invariant lines through the vertex of the pencil. In fact 
this hypercone is the only quadratic primal on which the surface 
can lie.* 

To determine the image of a plane field of lines in S, consider 
the S; which is determined in Sy by 7 and the plane, p, which 


* Cf. Baker, Principles of Geometry, vol. II, pp. 275, 279. 
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represents the lines of the given plane field. This S; intersects 
V,.° in the plane p, and in a V+‘, the Veronese surface, which is 
the transform of p. Now consider the Ss whose intersection with 

’ represents the special linear complex of lines meeting any 
one of the planes determined in S, by the points where an S; 
through the plane field intersects Z;, Le, and L3. This Ss con- 
tains ];, 2, 1s, and the plane p. It therefore-contains the S; de- 
termined by 7 and ~, and hence contains V+‘. All points of V+ 
are thus representations of lines of S; which meet all the planes 
determined by the respective triads of points in which the S3’s 
through the plane field intersect Zi, Lz, and L3. As a point locus 
this system of lines is a V;* having a plane of double points.* 


7. Analytic Procedure. Let the three fundamental lines in S, 
be determined by the respective pairs of points: (10000) (01000), 
(00100)(00010), (00001)(10100). The common transversal of 
these lines is the line joining (10000) and (00100). The Grass- 
mann coordinates of these lines, as read from the matrices of the 
points which determine each line, are 


Pye Pis Pu Pis Pes Poy Po Psy Pos Pos 
0 0 0 0 0 0 0 0 0 
0 0 0 0 0 1 0 
| ee ae 0 0 1 0 0 0 0 1 0 
fee eee | 1 0 0 0 0 0 0 0 0 


The quadratic identities which define V,° are 


— + PosPs = 0, 
Py3P 45 — PuPss + PisPss = 0, 
— + PisPxs = 0, 
PywP3 — Pi3P2 + PisP23 = 0, 
— + PisP23 = 0. 


The S; determined in Sy by i, kb, 1, and a general point 
l=(he-- - of can be written parametrically: 


y = Ah) + we) + + 


This meets V,° in the points h, 2, i, 1, and a fifth point 1’, given 
by 


* C. Segre, Enclykopaidie der Mathematischen Wissenschaften, Band III, 
2, Heft 7, p. 952. 
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A = — — leslss), 
(A) — lesa) , 
56 = as , 
4s. 
Thus the equations of the transformation are: 
= 134, bed los) 45, 
kh = 45, = loglo? 14s, 
(B) = = , 
= — 14), Riss! = + 12s), 


From equations (A) it is evident that the complex of invariant 
lines is given by 
— = 0. 


From (A) it also follows that the image of LZ; is any line of the 
special complex /,,=0; the image of Lz is any line of the special 
complex /,,=0; and the image of LZ; is any line of the special 
complex /.,=0. 

The singular lines are lines of the following systems: 


1. Lines for which J, = 25 = 1s = 0, (lines meeting 7 ). 
2. Lines for which = = = 0, (lines meeting 
2.1. Lines for which hy, = = = 0, (lines meeting Ly»). 
2.2. Lines for which /;2 + le3 = les = his — 1s = 0, (lines meeting 


3. Lines for which jj, = hs = 135 = 145 = 0, (lines lying in the 
S; determined by 
L, and Lz). 

3.1. Lines for which /\4 = Jey = /34 = 145 = 0, (lines lying in the 
S3; determined by 
L, and L3). 

3.2. Lines for which /,2 = Jes = Jog = Jo = 0, (lines lying in the 
S; determined by 
Le and L3) 
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SOME SYMBOLIC IDENTITIES* 
BY L. I. NEIKIRK 


Differential equations were first solved by symbolic methods 
in England and on the continent in the first half of the last 
century. The differentiation symbol was treated as a symbol of 
quantity with restrictions. Then followed symbolic treatment of 
invariants and covariants, Cayley’s hyperdeterminant, and 
Aronhold’s symbolic notations. These were followed by Blis- 
sard’s umbral notation in the theory of numbers. 

This paper is devoted to showing that these are all reducible 
to symbolic differentiation. 

If we represent differentiation by the symbol D and separate 
the symbols of operation from the symbols of quantity, then 
any analytic identity, such as ,(y) = ®2(y), will give an opera- 
tional identity, ®,(D) =®,.(D). 

If this operational identity is applied to a second identity 


F, (x) = 


the result will be a new identity. Most identities obtained in 
this way are easily obtained otherwise. The following are some 
examples. 

Invariants and covariants. If D,=0/dx; and Dz=0/0x2, then 


n n—r 
DoD, (a,x; + aexe) = nla, ao, 


where a" is a special form of degree u, while the operation on 
the general form gives 


r_n—r n n—1 
f(x, x2) = (aoxi x2 +--+) = 
We now transform our coordinates: 


or 


1 1 
= — — mX2), Xo = — (— fox, + £142), 
A A 


* Presented to the Society, June 18, 1936. 
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where 
m1 


£2 
Then if and 8/0X2=Dk, we write 


A= ~ 0. 


Ox, OX, Ox, Ox, 


in symbolic form 
Dif = (Dine — Dok) f, 
also 


D.f = Dim + Dek2)f; 
and if D’=0/dy we have 
(DD’) = (D\Df — 
[(Din, — Dat.)(— Din + 
— (— Dim + 2 — Df £2) | 
| D, D{ | D, D{ . 


n2 — & 
D. AID, Di 


For example 


gives 

Di} 

while 

21D. Dz 
Therefore 


A oA, — 4A 1A3 + 3A 2” = A*(aoa4 = 4a\a3 3d2”). 
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This is Cayley’s hyperdeterminant notation. This would be, in 
the Aronhold symbolic notation, 
1 | D, Di 
Ag 


1 


Blissard’s umbral notation. Let D=d/dy. Then the symbolic 
form of MacLaurin’s theorem is 


D? 


y=0 


+ DF) 


P= 


F(x) = FO) | 


or 


2 


x?D* 
F(x) = (1 + «D+ +.--- | 
y=0 


Now if 
Bey* 
My 3! +-:---, (e 
then 
FO) | =1; DF) = | = 
y=0 y=0 y=0 
and 


n(n — 1) 
(1, + nB,+ Bs) — B,, 


= [(1 + B)" — B*] =0 


becomes 
— 1 


= {[(1 + D)" — = 0; 


y=0 


and one of the principal identities used in Blissard’s theory be- 
comes 


= 
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af(x) 
dx 


{ + (D — — f(x — OD) |F(y)} = 0 


Blissard’s remark, “An equation which has a representative 
quantity is not susceptible to any algebraic operation by which the 
indices would be affected,” becomes 


(Df)? 


UNIVERSITY OF WASHINGTON 


ON FOURTH ORDER SELF-ADJOINT 
DIFFERENCE SYSTEMS* 


BY V. V. LATSHAW 


A linear difference expression for which the differential trans- 
form is self-adjoint (anti-self-adjoint) we shall call self-adjoint 
(anti-self-adjoint).t We choose two fourth order difference equa- 
tions 


L*(u) = p(x) [u(x + 2) + u(x — 2)] 


1 
(1) + r[u(x + 1) + u(x — 1)] + R(x)u(x) = 0, 
L-(u) = p(x) [u(x + 2) — u(x — 2)] 

(2) + r[u(x + 1) — u(x — 1)] = 0, 
where Lt+(u) is self-adjoint and L~(u) anti-self-adjoint for the 
range (x=a,a+1,---,b—1;b—a2=4). R(x) and p(x) are both 
real, p(x) being a non-vanishing periodic function of period two; 
d is a parameter. 

Let the functions (1, ye, v3, ys) constitute a fundamental set 
of solutions for either (1) or (2), and (wi, we, w3, wy) the set ad- 
joint to it. The two sets are related by the equations 


* Presented to the Society, October 30, 1937. 

+ J. Kaucky, Sur les équations aux différences finies qui sont identiques a 
leurs adjointes, Publications of the Faculty of Sciences, University of Masaryk, 
No. 22 (1922). For a discussion of adjoint differential expressions of infinite 
order, see H. T. Davis, The Theory of Linear Operators, 1936, pp. 474-475. 


= 
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+ 1) 
w;(x)p(x) = 
(3) yi(x — 2)--- ya(x — 2) 
A(x) = ; 


+ 1) -- + 1) 
in which A (x) denotes the Casorati determinant and A ;;(x) the 


cofactors of its elements. Since w; satisfies the equation L(u) =0, 
one can write 


(4) wi = cayi(x) +--+ + cisya(x), (i = 1, 2, 3, 4). 
The c;; have properties stated in the following theorem.* 


THEOREM 1. The matrix C* of the substitution (4) is skew-sym- 
metric and C~ symmetric. Further, there is a set of relations involv- 
ing c;* and the second order minors of A*(x) from which the c;* 
may be calculated explicitly: 


4 Vi Vi 1 


1) ~ p(x)p(x + 1) 


4 Vi Vi 1 
1 ‘) 1 1 p(x + 1) 


We resolve (4) to get 
1 
(5S) [wi(x — + --- + whe + 1)A,,(x) | 
A(x) 


We use (1, 2, 3) to furnish the following relations 


* We introduce a (+) and (—) convention to distinguish between quanti- 
ties associated with (1) and (2) respectively. Also we set 


o\t 
u(x + r)o(x + s) + v(x + r)u(x + s) = 
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A4;(x) 
w(x — 1)p(x — 1) = i? 
w(x — 2)p(x) + — 1) = 
A1;*(x) 
+ A2;*(x) 
+ p(x + + = 


Combining (5) and (6) we find 


+ wt w*\t wt 
= ple + A ) + ) 


q 1 2 
(7) 
0 --1 
Obviously we have c,;j+= —c;# and cj-=c;7. Combine (3) and 


(4) to eliminate the w’s. These equations together with L+(y+) 
=0 yield the set of equations involving the minors. The sixth 
order determinant composed of the two-rowed minors is non- 
vanishing (= [A+(x) ]*).* 

In the development of adjoint difference systems the formula 


(8) f(x + 1) = f(x) + f(x) 


is used to provide a Lagrange relation 
b—1 
(9) [oL(u) — uL(v)] = = +--- + UsVi, 


in which the U; are eight linearly independent forms arbitrarily 
chosen: 


(10) U; = ayu(a — 2) + + 1) 
+ byu(b — 2) + + 1). 
If the two systems 
L(u) = 0, U, = U2 = U; = Ui = 0,7 


(11) 
=0, =. =0, 


* Turnbull, The Theory of Determinants, Matrices, and Invariants, 1929, 
p. 87. 


z=a 
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are equivalent, we shall call them self-adjoint difference sys- 
tems. 


THEOREM 2. Given the fourth order difference systems composed 
of the equations 
L+(u) = 0, U; =0, (i = 1, 2, 3, 4), 


defined in (1), (2), (10), let u(x) and v(x) be any pair of functions 
satisfying U;=0; then II(u, v)=0 is a necessary and sufficient 
condition that the given system be self-adjoint. 

That the above condition is necessary needs no proof. Let 
(1, U2, us, us) be four linearly independent functions satisfying 
U;=0 and u any linear combination of them. Through substitu- 
tion the identity II(u, v) =0 gives 


(12) Us(u)Va(u) + --- + Us(u)Vi(u) = 0, = 1, 2, 3, 4). 


Since the set (Ui, U2, ---, Us) is linearly independent, it fol-- 
lows that the four systems of constants comprising the coefh- 
cients of the V’s are linearly independent and we have 
Vi(u) = Ve(u) = V3(u) = Vs(u) =0. A similar argument shows 
that any function satisfying the given boundary conditions 
will also satisfy the adjoint boundary conditions. 
We record some examples which fulfill the condition for self- 

adjointness. 

= u(a — 2) + aiu(a — 1) + p(a)u(a) 

bi3(b) bi4u(b 1) = 0, 
Ust = + ayou(a — 2) + u(a — 1) + du(a) 

+ p(a + 1)u(a + 1) + bogu(b + 1) = O, 
Us = + bizsu(a — 2) + p(b)u(b — 2) + Au(b — 1) 

+ bszu(b) + bsu(b + 1) = 0, 
+ — 2) + boyu(a — 1) + p(b + 1)u(b — 1) 
+ bsyu(b) + dif u(b + 1) = 0, 


+ I+ 


with the agreement that =0. For self-ad- 
joint Sturmian boundary conditions we make the added restric- 
tions, bis = bis = bos =(. 

We now introduce a function G+(x, t) defined for (a—1<x<b), 
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(ast<b), satisfying the given boundary conditions and for 
which L+[G+(i, j) ] =6;;. 

THEOREM 3. Let dj*,d2o*, - - - be sets of characteristic values for 
the systems L+(u) =0, U;=0 defined in (13). There exists a dit 
in the interval (\7 SAF provided 

G-(b, b — 1) —Gt(b, b — 1) — [G-(« — 1, x) + Gt(x —1, x)] #0, 
(r Sd SdH). 


If (6—a) ts an odd integer, then every value of \~ is a characteristic 
value. 


Let D+(\)=0 be the characteristic equations for (13). By 
writing D+(\) in determinant form one finds 


(14) — = 2D+ | b-—1)+ 


z=a 


This relation enables us to write 


d {D~ 2D- 
xe 
dy \D* Dt 


(15) 
— — 1, x) + G(x 1, 

Between the two real consecutive zeros \7 and d;,, of D-(A) 
either the bracketed expression or D+(A) must vanish. By as- 
sumption the bracketed expression does not vanish. 

The proof of the final statement in the theorem consists 
merely in noticing that D~(A) for this case is a skew-symmetric 
determinant of odd order. 


LEHIGH UNIVERSITY 
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W. E. ROTH 


ON CERTAIN MATRICES AND THEIR 
DETERMINANTS* 


BY W. E. ROTH 


1. Introduction. It is well known that certain determinants, 


| cs; , (i, 7=1,2,---, 1), are expressible as the products of r 
linear factors of the form adjwatapwet --- +4;,0;, (4=1, 
2,---,7), where a;; are rational functions of the elements, ¢;;, 


in the first row of |c;;|, and where w,;; are functions of certain 
roots of unity and depend only on the relations existing between 
the elements of succeeding rows and those of the first row of the 
given determinants. Cyclic determinants and some related types 
which have this property are well known.f 

In the present paper a generalization of these results is carried 
forward in two directions; first, we replace the scalar elements 
in the array of the given determinant by square matrices of 
order n; and second, we permit succeeding rows of the array to 
be other than permutations of the matrices occurring say in the 
first row of the given determinant. Certain special types, which 
we here generalize, have been studied by Puchta, Noether, 
Baltzer, Drude, Burnside, Scorza, and others.{ 


* Presented to the Society, April 9, 1937. 

+ Pascal, Die Determinanten, §§20-21. Gegenbauer, Uber eine specielle 
symmetrische Determinante, Sitzungsberichte, Akademie der Wissenschaften, 
Mathematische-Naturwissenschaftliche Klasse, vol. 82, II-III (1880), pp. 
938-942. Burnside, On a property of certain determinants, Messenger of Mathe- 
matics, vol. 23 (1894), pp. 112-114. 

t Puchta, Ein Determinantensatz und seine Umkehrung, Denkschriften der 
Wiener Akademie, vol. 38, 2*e Abth. (1878), pp. 215-221; Ein neuer Satz aus 
der Theorie der Determinanten, ibid., vol. 44, 2'e Abth. (1882), pp. 227-282. 

Noether, Zur Theorie der Thetafunctionen, Mathematische Annalen, vol. 16 
(1880), pp. 322-325; Notiz iiber eine Classe symmetrischer Determinanten, ibid., 
pp. 551-555. 

Baltzer, Ueber einen Satz aus der Determinantentheorie, Nachrichten, 
KG6niglich Gesellschaft der Wissenschaften, Gottingen, 1887, pp. 389-391. 

Drude, Ein Satz aus der Determinantentheorie, ibid., pp. 118-122. 

Burnside, loc. cit. 

Scorza, Sopra una certa classe di determinanti e sulla forma Hermitiane, 
Giornale di Matematica, vol. 51 (1913), pp. 335-342, 
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2. Direct Product Matrices. lf A and B are nXn matrices and 
if 
U= V= (v4), (i,j = r), 


are r Xr matrices whose elements u;; and v;; are commutative 
with those of A and of B, then the mr Xur matrix 


(A)U U(A) (u;;A) = (Au;;), 
is the direct product of A and U, and it readily follows that 
(1) (A)U-(B)V = (AB)(UV) = (UV){AB). 


Moreover, the matrices (A )U and A(U) are equivalent in that 
their characteristic matrices have identical elementary divisors.* 


3. The Matrix (Ao)I,+(A1)U+ --- +(A,1)U™. We shall 
assume that the elements of the Xn matrices, A;, (¢=0, 
1,---,r-—1), belong to the field F, and that the r Xr matrix U, 
likewise with elements in F, is such that U—XI has the elemen- 
tary divisors (t=1,2,- - -,s), where (@=1,2,--- ,7r), 
are not necessarily in F. Let the non-singular matrix, Q, with 
elements in F extended by the adjunction of the roots of 
| U—XI| =0, be such that 


(2) QUO = Ut 
where 
i 0- 
6; 1 
(3) U;= 
0 0 


is an 7;Xr; matrix and f=. For convenience, we shall 
designate as B(U) the ur Xnr matrix discussed in this section: 


where J, is the unit matrix of order r. Hence according to (1) 
and (2) we have 


* For further details on the algebra of direct product matrices see MacDuf- 
fee, The Theory of Matrices, Springer, 1933, Chapter VII; and Roth, On direct 
product matrices, this Bulletin, vol. 40 (1934), pp. 461-468. 


— 
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= B(U:) + B(U:) + BID, 


and 
BOG) B’(6;)--- ) 
| B(@;) --- 
oO --- B® J 
where 


(k=0, 1,---,7r;—1; 2,---, 5), the C’s being binomial 
coefficients. We have proved the theorem: 


THEOREM. I[f 
B(U) = (Ao)I, + (A1)U +--+ + 
where the A; are nXn matrices with elements in F, and if U, an 
rXr matrix with elements in F, is equivalent to the direct sum 
U,+U2.+ --- +U.,, where the U;, (i=1, 2, --- , 5), are given by 
(3); then B(U) ts equivalent to the direct sum 
B(U,) + + --- + B(U,), 

where the B(U,), (i=1, 2,---, 5), given by (4) have their ele- 
ments in F extended by the adjunction of the roots of | U—nxI| =0. 

This theorem, whose proof as given above is almost trivial, 


has very extensive applications in the theory of determinants. A 
few of these will be illustrated in the following section. 


4. Determinants. The theorem above leads at once to the fol- 
lowing: 

Coro.vary. If (i=1, 2,--- , 1), are the characteristic roots 
of U, whether distinct or not, then 


(5) | B(U)| = II | | .* 


i=l 


* Williamson (The latent roots of a matrix of special type, this Bulletin, 
vol. 37 (1931), pp. 585-590, Theorem I) and Finan (A theorem on matrices, 
abstract 41-11-359, this Bulletin) have obtained results closely allied to this 
corollary. 


= 
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That is, the mr Xur determinant | B(U)| is expressible as the 
product of r determinants, | B(O,)|, (i=1,2,---,7), of order n. 
For example, the determinants studied by Baltzer* and by 
Scorzaf would, according to our notation, be given by 


A | Ag A; | | ‘Ao)I h U ( 0 7 
and since | U—\J| =(i—d)(—i—d), we have 

A =| Ao + iAy| -| Ao — |. 


Hence the given determinant is the sum of two squares provided 
the elements of Ao and A; are real. Also we can conclude that if 
Ao+tA;, where Ao and A; are real, have a real characteristic 
root, then (A»9)J2+(A1)U above has this root as a multiple char- 
acteristic root. 

The cyclic determinant 


| 
where 
(oO 1 0---0) 
10 0 1---90| 
lo 0 0---1 
0 
and the A,;, are matrices, is given by 


the product 


II 
i=1 


* Baltzer, loc. cit. 
+ Scorza, loc. cit. 
t Wedderburn, Lectures on Matrices, p. 101. 


— 
— 
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where w is a primitive rth root of unity, that is, a primitive root 
of | U—XI,| =0. 

Each of the factors, | B(6,)|, (i=1, 2,---, 7), in (5) is again 
expressible as the product of s factors, if the Xn matrices Aj, 
(i=0, 1,---, r—1), are given by 


A; = + (Aia)V + + 


where the A ,;, (¢=0, 1, -- - ,r—1;7=0,1,---,s—1),aremxXm 
matrices, where V is an sXs matrix, and »=ms. This factoriza- 
tion may be continued still farther under certain conditions 
which are now obvious. 
Matrices of the form given by B(U) above have arisen fre- 
quently in the literature, particularly in the theory of algebraic 
numbers. For example, the “Xn matrix whose elements are 
algebraic numbers with 9, a root of the equation 
Go tary+--- =0, 

as basis, can be written 
B(@) - Ao + A 10 + gt + Ae, 

where the A;, (i=0, 1, - - - ,r—1), are Xn matrices. Then 
6B(@) = Ad A {6 eee 

where 


Similarly 
= Ad’ + + 4,107", 
where 
Ad’ Af’ (i= 1,2,---,7r—1), 


and so on. Hence the ur Xn matrix 


BQ) | 1 


| (r—1) (r—1) (r—1) 


| 
} 


E 
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0 
= BCU): {Ir)| |, 
where 
bs 0 
| 0 
(6) Ya 
0 0 0 1 
do a4) 
and 


B(U) = (Ao)I, + (A1)U +--+ + 


The ur Xnr matrix B(U) in F corresponds to the »Xn matrix 
B(6), whose elements are algebraic numbers, in the sense of an 
isomorphism under addition and multiplication.* 


UNIVERSITY OF WISCONSIN, 
EXTENSION Division, MILWAUKEE 


* The determinant, B(U), may be regarded as the norm of B(@). See 
Dickson, Algebras and their Arithmetics, p. 70. 
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REMARKS ON THE INITIAL VALUE PROBLEM OF 
THE GENERAL PARTIAL DIFFERENTIAL 
EQUATION OF THE FIRST ORDER 


BY H. R. COOLEY 


1. Introduction. This paper presents a new approach to the 
classical theory of partial differential equations of the first order 
by establishing the equivalence of the initial value problem to 
another problem associated with a system of quasi-linear partial 
differential equations of the first order. All the equations of this 
quasi-linear system will have the same “principal part”; this 
means that the derivatives have the same coefficients in all the 
equations. Incidentally, systems of differential equations with 
the same principal part are assuming an important réle in the 
theory of partial differential equations.* 

In the case of quasi-linear systems the theory of integration 
can, as we shall see, be reduced to that of the integration of a 
system of ordinary differential equations in the same way as 
with a single linear partial differential equation. We shall then 
show that the initial value problem of the general partial differ- 
ential equation of the first order is equivalent to a certain initial 
value problem of our quasi-linear system with identical princi- 
pal parts. In this way we can develop the theory of characteris- 
tics and the solution of the initial value problem for a general 
partial differential equation of the first order. 


2. Systems of Quasi-Linear Partial Differential Equations of 
the First Order with Identical Principal Parts. Let 


(1) a;— = (u = ,m), 


be a system of quasi-linear partial differential equations in 
which the coefficients a; are the same in all equations. This 
is a system of quasi-linear equations with identical principal 
parts for m functions ™, , Um of k variables x, --- , x,. In 
these equations a; and b, are functions of x1, - - - , X%, M1, °° Um. 


* For this topic see Courant-Hilbert, Methoden der Mathematischen Physik, 
vol. II, Chapter 2, appendix and Chapter 5 (in press). 


862 
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If now, generalizing a well-known procedure,* instead of re- 
quiring m functions u, explicitly, we require m implicit relations, 


$.(%1, Um) = Cu, (u = 1, m) , 


depending on m parameters c,, the system (1) is transformed 
into a system of homogeneous linear partial differential equa- 
tions with identical principal parts, or, what amounts to the 
same thing, a single homogeneous linear partial differential 
equation for which we require m independent solutions. For 
these functions ¢,, we have 

0d, ™ 


(2 — + =0 
) Ox; Ox; 


Multiplying by a; and adding, we obtain 
2 0%, 


i=1 Ox; Olly Ox; 


(3) Dd db = 0. 


Ox; OM 


If we let =O, R-+-m =r, equation (3) be- 
comes 
r 0d 
(4) a, — = 0. 
OX, 

Equation (4) is a homogeneous linear partial differential equa- 
tion and is equivalent to the system of ordinary differential 
equations 

dx, do 


= 


0, 
ds ds 


which are its characteristic equations. By integrating this sys- 
tem of characteristic equations a solution of (4), depending on r 
independent parameters, can be found. 

In order to show the equivalence of system (1) and equation 
(4), in addition to the fact that a solution of (1) satisfies (3) in 


* Cf. Courant-Hilbert, loc. cit., Chapter 1. 


Since ,a;0u,/dx;=by, we have 
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accord with the above derivation, it is necessary to show that, 
conversely, m independent solutions of (4) furnish a solution of 
(1). 
Let ¢1,---, 6m be m independent solutions of (4). Then 
= |d¢, /Auy| ~0, and hence from the equations ¢, =c, we can 
obtain 


These functions u, must satisfy (1). For, differentiating 


we obtain (2). Using this equation with (3), we have 


A=1 i=1 Ox; O86 ax, 
or 
0g = Ou 
Hence 
m do 
h=1 Ox; 


Since A0, it follows that d) >. ,4,0u,/dx;=0, that is, (1) is 
satisfied. | 


3. The Derivation of a System of the Form of (1) from the Gen- 
eral Partial Differential Equation of the First Order. If the num- 
ber of independent variables is k, the general partial differential 
equation of the first order is 


: Ou 
(5) Xk, u, » Pr) = 0; 2; = ——* 

Ox; 
Differentiating (5) with respect to x,, and making use of the 
relations 0p;/0x;=0p;/0x;, we obtain the first k equations of the 
system 


f k OF OF OF 
a — — = — —-— py, =1,---,&), 
() Op; Ox; Oxy Ou 


(6) 


k OF ou k OF 
b = — p;. 
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The last equation is a formal triviality depending on the rela- 
tion 0u/dx;=p;, written down to complete the system. This sys- 
tem is of the form (1), in which 


x; = Xi, Uy, = Py, ) 
OF OF oF (u=1,---,k). 
Op; Ox, Ou 
= u, => pi- 
int OP; 
The characteristic equations in this case are 
dx; OF du > OF OUy OF OF 
ds Op: = Op: ds Ox, Ou 


4. Solution of the Initial Value Problem of the General Partial 
Differential Equation of the First Order. The initial value prob- 
lem for equation (5) is the following: We are given an initial 


manifold, M(x1, x2,---, x.)=0, with a set of functions u, 
fi, Pe, satisfying the relations 

(7) du = >> pidxi, 

i=1 


which, for 0\//dx;,0, is equivalent to the k—1 relations 
(— ) (— ) <b) 


(8) F(x, pi; px) = 0 


on M=0. These initial data are assumed to be non-characteris- 
tic; this means that for the initial values 


OF 0M 
(9) 
Op; Ox; 
holds. Then we seek a function u(x:,--- , xx) satisfying the 


differential equation (5), which has the values prescribed for u 
on M, and for which the derivatives 0u/0x; have the prescribed 
initial values p; on M, 


(i,4=1,---,). 

and 
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Again in the initial value problem for the system (6) the val- 
ues of u and p; are prescribed on the initial manifold M, subject 
to the same condition (9) in order that M be non-characteristic. 

Now our statement is as follows: If in this second problem the 
additional conditions (7) and (8) are satisfied, then this problem 
becomes equivalent to the first problem for equation (5), and by 
solving this second problem for the quasi-linear system we have 
solved the problem for the general equation. 

The equivalence of the initial value problem of the general 
equation (5) and the special initial value problem of the system 
(6) can be shown in the following manner. 

First we note that a solution of the initial value problem of 
(5) yields a solution u, p;=0u/0x;, of the corresponding initial 
value problem of the system (6), since if u=u(x;, - -- , x.) satis- 
fies (5) identically, then (6) follows immediately by differentia- 
tion. 

Now let u(x, --- , xx), , Xe), in whichi=1,---, k, 
be a solution of the restricted initial value problem of (6). We 
first observe that from (6b) and (7) we can infer that on the 
initial manifold M the relations 0u/dx;=p; hold, since the de- 
terminant of the linear equations (6b) and (7) for 0u/dx;— p; is 
(except for a factor) the left side of (9), and thus does not van- 
ish. Furthermore from (8) we have F=0 on M. We now have to 
show that F=0, 0u/dx;= pi, everywhere. 

Let 


Ou 
ax; Pi = P(x, Xx). 


We have R=P;=0 on the initial manifold. Differentiating R 
with respect to x;, and rewriting equations (6), we have 


(OR OF ou OF 
Ox; 


Ox; Ox; 


OF k OF Op; 
(10) 0=—+4+F.p:+ — (@=1,---,hk). 
Ox; ju OP; O%; 
k OF ou 
| 
OP; j=1 OPj OX; 
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Multiplying these 2k+1 equations respectively by 
OF OF OF OF 
Opi Op: Op, Opi. 


F,, 


and adding, we obtain 


This is a homogeneous linear partial differential equation for R, 
with R=0 on the initial manifold. Therefore R=0. 
The last of equations (10) can be written 

(11) > — P;=0. 

j=1 OP; 
Making use of the fact that 0R/dx;=0, and subtracting the 
(k+7)th equation from the ith equation, in system (10), we 
obtain the system 


0 k OF (0p; Op; 


ax; j=1 OP; \Ox; Ox; 
or 
k OF fOP; oP; 
jut OP; \O%; Ox; 


If we differentiate (11) with respect to each of the independ- 
ent variables and substitute in (12), system (12) then becomes 


OF oP; 
(13) =A, 
j=1 OP; OX; 
in which A; is a homogeneous linear form in Pi,--- , P,. The 


left-hand side of each equation of (13) is the directional deriva- 
tive of P;, 0P;/dc, in the direction determined by 


where o is a parameter of the characteristic curve through the 
point under consideration. Hence (13) may be written 


k OF AR 
> — — -0. 
im OP; O2; 
Ox; OF 
be 
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(14) 

This system, together with the initial conditions, is satisfied by 
P;=0, (¢=1,---, Rk). Hence, on account of the uniqueness of 
the solution of (14) with given initial values, we conclude that 
P;=0, and the proof is complete. 


New York UNIVERSITY 


ON THE EXISTENCE OF LINEAR FUNCTIONALS 
DEFINED OVER LINEAR SPACES* 


BY R. P. AGNEW 


1. Introduction. A function g(x) with domain in a linear space 
E and range in the set R of real numbers is called a functional, 
and q(x) is called linear, if 
(1) + by) = ag(x) + bg(y), yeE; a,beR. 
We call a functional r(x) an r-function (over £) if there exists a 
linear functional f(x) with 


(2) f(x) S r(x), xeE. 
Using a notation of Banachf we call a functional p(x) a p-func- 
tion if 

(3) p(tx) = tp(x), t= 0, xeE, 
(4) p(x + y) S p(x) + ply), «,yeE. 


A fundamental theorem of Banach (loc. cit., p. 29) can be 
stated as follows: 


THEOREM (Banach). Each p-function is an r-function. 


In some problems{ involving existence of linear functionals 
fi(x) having prescribed properties, there appears a functional 
q(x) with the following significance: There exists a linear func- 
tional f; having the requisite properties if and only if there exists 


* Presented to the Society, September 8, 1937. 
7 S. Banach, Théorie des Opérations Linéaires, Warsaw, 1932, p. 28. 
} The author intends to discuss these problems at some future time. 


A 
| 
| 


1937-1 LINEAR FUNCTIONALS 869 


a linear functional f with f(x) <q(x), that is, if and only if g(x) 
is an r-function. If g(x) can be shown to be a p-function, the 
problem is solved by Banach’s theorem; if g(x) is not a p-func- 
tion or one is unable to decide whether g(x) is a p-function, 
Banach’s theorem cannot be applied. These considerations, and 
the fact that it is easy to give examples of r-functions which are 
not g-functions, lead one to desire an analytic characterization 
of r-functions. In §2 we give such a characterization, and in 
§3 we give some closely related theorems. 


2. Characterization of r-functions. We prove now the theorem: 


THEOREM 1. In order that a functional r(x) defined over E 
may be an r-function, it is necessary and sufficient that 


>0 


(5) ‘ie 


ty 


In (5), >ox, stands for the sum x;+ --- +x, of elements 
x,eE. To prove necessity, let r(x) be an r-function and let 
f(x) be a linear functional with f(x) <r(x) for all xeZ. Then if 


n, h, h, >0 and > x. =0, we have 
(6) (xe) = f(texe)/te S r(texe)/te, 
so that 


(7) 0=f0) w = 


and (5) follows. 
To prove sufficiency, let (5) hold and define the functional 
p” (x) by the formula 


(8) p(x) = xeE. 


ty 


To show that p‘” (x) exists (is finite) for each xeE, we observe 
that if m, 4,---, and }ox,=x, then m+ --- +x, 
+(—x) =0 and it follows from (5) that 


0 


IV 


> r(— x) 
k=1 ty 1 
and hence 


R. P. AGNEW 


k=1 ti 
which implies that —r(—x) Sp” (x). If in the sum in the right 
member of (8) we put =1, 4: =1, x1 =x, we obtain (x) <r(x). 
Therefore 


(9) — r(— x) S p(x) S r(x), xeE. 


We prove next that p‘” (x) is a p-function. If xeE and ¢>0, 
then 
= 


cis. 


kml ty 


gis. 


pr (tx) 


n th, 
ea, 
k=1 


so that (tx) =tp' (x) for t>0. Substitution of x=0 
in this formula gives p‘” (0) =0. Therefore 


(10) = (x), t= 0; 
To prove that 
(11) p(x + y) p(x) + p(y), x,yeE, 


let x,yeE be fixed and let e=0. Choose m, hh,--- , tm>O and 
XmeE such that }>x;=x and 


m 


j=l 


and choose m, m,---, and yi,---, y,e£ such that 
> y.=y and 


n 


< p(y) +. 

k=1 
Since m+n, t;, uz, >O and 4+ --- +xmtyit: tyn=xt+y, 
it follows from the definition of p (x+y) that 


870 (December, | 
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j=1 i k=1 Uk 

The arbitrariness of gives (11). Thus (x) is a £-function 

and it follows from Banach’s theorem that there exists a linear 

functional f(x) with f(x) (x). Using (9), we obtain 

f(x) Sr(x); thus r(x) is an r-function and Theorem 1 is proved. 


3. Significance of p‘” (x). From Theorem 1 and its proof, 
we obtain the first part of our next theorem. 


THEOREM 2. Jf r(x) is an r-function, then the functional 
p” (x) defined by 


is a p-function with 

(13) — r(— x) S — p(x) p(x) r(x), xeE; 
moreover if p(x) is a p-function with p(x) <r(x) for all xeE, then 
(14) — x) S — p(— S p(x) pM(a), — we. 


In establishing (13), we use (9) and the fact that, for any 
p-function, 0 = p(0) = p(x —x) p(x) + p(—x) and hence — p(—x) 
<p(x) for all xeE. If p(x) Sr(x); n, t,---, t,>0; and 

x,=x; then 


p(x) = p( xx) = pltexr)/te < /ty 


and p(x) <p” (x) follows. The remaining inequalities in (14) 
follow easily, and Theorem 2 is proved. The gist of Theorem 2 
is that p( (x) is the “greatest” p-function p(x) with 
p(x) Sr(x). In particular, if r(x) is a p-function, then p‘” (x) 
=r(x). 

Since each linear functional f(x) is a p-function, Theorem 2 
implies the following theorem: 


THEOREM 3. If r(x) is an r-function and f(x) ts a linear func- 
tional with f(x) <r(x), then 


(15) 2) S — 2) S f(x) S p(x) S x(x), 


n n n 
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It thus appears that the class of linear functionals f(x) for 
which f(x) Sp” (x) is identical with the class of linear func- 
tionals f(x) for which f(x) Sr(x). 

4. Conclusion. The functionals g(x) mentioned in the intro- 
duction often have the property q(tx) =tq(x) for 120, and xeL£. 
Hence it is of interest to note that if 
(16) r(tx) = ir(x), +20, xeE, 
then the criterion (5) that r(x) be an r-function reduces to 
(17) g-l.b. r(x.) = 0, 


and that formula (12) for p‘” (x) reduces to 


n 


(18) = g.l.b. > r(x), 


k=1 


CORNELL UNIVERSITY 
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ON REFLECTION OF SINGULARITIES OF HARMONIC 
FUNCTIONS CORRESPONDING TO THE 
BOUNDARY CONDITION 0u/dn+au=0 


BY HILLEL PORITSKY 


1. Introduction. Familiar “reflection” principles across a plane 
at which a harmonic function u satisfies either of the two bound- 
ary conditions 


(1) u=0, 
(2) du/dn = 0, 


where 0/0n denotes the normal derivative, extend the function 
from one side of the plane to the other one by means of its nega- 
tive or positive image respectively. In particular, the singulari- 
ties of u to one side of the plane are also reflected into their 
negative or positive images. 

In the following we consider the nature of the “reflection” or 
analytic continuation of a harmonic function u across a plane 
boundary corresponding to what is perhaps the next simplest 
boundary condition, namely: 


(3) —+au=0, 


where a is a constant. It is shown that the image of each singu- 
larity So of u is relatively complex and consists of 

(a) a positive image S, of So in the boundary plane; 

(b) an exponential trail of negative images along the line 
through Sp and S;, beyond S,, and totalling in amount 
double the negative of So. 

Results similar to the above are established for other differential 
equations; for instance, for the equation of heat conduction. 
Conditions with higher order derivatives are also considered. 
Aside from the interest of the subject matter in connection 
with analytic continuation, as well as from the point of view of 
general curiosity that makes one “peep behind the looking 
glass,” the subject is also of interest in view of several pos- 


| 
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sible applications. One application is to the problem of heat 
flow from an underground source; another one is to the problem 
of eddy currents induced in a semi-infinite solid with a plane 
boundary by nearby alternating currents. These applications 
are discussed. 


2. Analytic Continuation of Harmonic Functions Correspond- 
ing to (3). Consider a function u, harmonic for x20 and satisfy- 
ing the boundary condition: 


Ou 
(4) —-—au=0 along x=0. 
Ox 
Let 
(5) ou 
w=—— an. 
Ox 


The function w is also harmonic, while along x=0 it vanishes. 
Hence w may be continued analytically to «<0 by means of 
negative reflection : 


(6) w(— x,y,z) = — w(x, y, 2). 


Substitution from (5) converts (6) into a differential equation 
whose solution yields 


z 


(7) u(— x) = — u(x) + 2a f dx’. 

0 
This is what the simple negative and positive reflection corre- 
sponding to (1) and (2) respectively is to be replaced by in case 
of (4). 


3. Green’s Function Corresponding to (4). To examine the na- 
ture of reflection of the singularities of harmonic functions satis- 
fying (4), consider for x >0 a function u which satisfies the fol- 
lowing requirements: 


((a) wis harmonic except near 0, 0). 
'(b) near Po, wis of the form u=1/ro+u’, where ro is the 
(8) 4 distance from Po, and u’ is harmonic. 
\(c) u and its first derivatives vanish at infinity. 
\(d) along x =0, u satisfies the boundary condition (4). 
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This function u may be designated briefly as the “Green’s func- 
tion” for x >0 with pole at Po, corresponding to the boundary 
condition (4). 

If the real part of @ is non-negative, 


(9) R(a) 2 0, 


a solution of (8) is given by 


(10) v = 1/ro + 1/r, — 2al, 
where 

—h 
(11) I -f (et dx’: 


here 1; is the distance to P,=(—hA, 0, 0) and r’ the distance to 
the point (x’, 0, 0) on the x-axis. If we use an electrostatic ter- 
minology which designates 1/ro as the potential of a unit charge 
at Po, the solution v may be described as the potential of unit 
charges at Po and P, and of a line distribution of charge of 
density — 2a e***+” per unit length of the x-axis, extending from 
x=—-« to P,. The total charge of the line distribution is 
—2a dx = —2. Added to the unit charges at Po and P; 
this results in a net charge zero. Thus v vanishes at infinity and 
to a higher order than 1/r. 

To prove that v satisfies the boundary condition (4), one 
proves 


1/r, + al 


by differentiating under the integral sign, replacing 0(r~') /dx 
by —0(r-')/dx’, and integrating by parts. It follows that 


Ov 1 1 0/1 1 
(12) = a(—-— +—|{—+—}. 
Ox ro Ox 


Since the latter vanishes at x =0 the proof is complete. 


4. Ways of Arriving at the Preceding Green’s Function. Its 
Uniqueness. Introduce the function w given by (5). The latter 
is harmonic for x >0 except at Po where it becomes infinite like 
0(1/ro) /8x—a/ro; it vanishes at infinity and also vanishes at 
x=0. It follows then that wu satisfies (12): 


— 
= 
_ 
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ou 1 1 
Ox Ox Yo Ox 


Integration of (13) in the form 
(14) «= -f y, z)dx’ -f e-**w(x + 5, y, 2)ds, 
z 0 


and consequent integration by parts and simplification, lead to 
(10), (11). 

Another way of obtaining u—and this is the way in which it 
was first obtained—is by means of the Fourier or Fourier-Bes- 
sel integral. This will be found useful for the case in which a 
does not satisfy (9). In the latter case it will be seen that while 
(13) still persists, the integrals in (11), (14) are divergent. 

Write u in the form 


1 
(15) “u=—-+ f dr, 
To 0 
where p?=y?+2?; utilizing the familiar ‘ategral expansion 
(16) 1/r = (p? + x?)-1/2 -{ J 
0 
and letting the resultant integrand of u satisfy (4), one obtains 
(17) f) dh Ah 2a 
=¢° 
A+ a 


The first term in (17) leads to the positive image at P;. The sec- 
ond term yields 


J Xp dr 
o(Ap)dr 
A+a 


By deforming the path of integration in the complex \-plane the 
analytic continuation of u to x <0 may be obtained and a charge 
density — 2a e*‘*+” deduced along the negative x-axis for x << —h. 
For negative a the path of integration in (18) must avoid the 
pole \= —a and hence cannot be confined to the real axis. We 
shall not stop over the particulars. The analytic continuation of 
harmonic functions by deforming the path of integration of 
proper integrals will be considered at a future date. 
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It is worth noting that for negative a the solution of (8) is no 
longer unique, since many harmonic solutions of the homogene- 
ous equation 0u/dx—au =0 (see (13)) exist which are harmonic 
and vanish at infinity for x >0, namely, 


(19) C(y, z)e*7, 

where C is a solution of 

(20) 
dy? dz? 


Thus, when we confine ourselves to axial symmetry, u is arbi- 
trary to within a multiple of 


To(ap)e**, 


which will be recognized as proportional to the residue of the 
integrand (18) at the pole X= —a. 

For R(a)20 the Green’s function vanishes at infinity for 
x=0; for R(a) <0, if the vanishing at infinity (condition (8d)) 
is to be interpreted as applying for all x, then mo solution of (8) 
exists. 


5. The Reflection of Singularities of Other Harmonic Functions. 
The singularities of other single-valued harmonic functions can 
be obtained by reflections similar to that of the Green’s func- 
tion. For instance, the harmonic function which becomes singu- 
lar at Po like 0(1/r0)/0y (while satisfying (8a), (8b), (8d)) is 
equal to 0v/dy, where v is given by (10). Describing the singu- 
larity at Po as a dipole of doublet of unit moment whose axis is 
in the direction of the y-axis, we see that the analytic continua- 
tion to x <0 has a similar dipole at P; and an exponential trail 
of dipoles along the x-axis beyond P;. This follows from the 
fact that if u satisfies (4), so does also 0u/0y, or it may be proved 
by replacing the derivative by a limit of a difference quotient, 
interpreting each of the terms of the difference as a potential of 
a point charge, applying the now familiar reflection to it, and 
passing to the limit. 

The harmonic function u which becomes singular like 
0(1/ro)/Ox at Po is not given by 0v/dx where v is given by 
(10) since, if « satisfies (4), du/dx will not in general satisfy it; 
this harmonic function may, however, be obtained by the limit- 
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ing process which analyzes the singularity at Po into singulari- 
ties of point charges. There results the formula 


1 1 
To 2a —h ea(zth) 
(22) u= — — + 2a’ da’, 


Ox Ox Lal 


, 
r 


exhibiting for u a dipole at P; but of opposite moment to that 
at Po, a point charge of amount —2a, and an exponential trail 
of positive charges. 

Harmonic functions which satisfy (4) but fail to satisfy 
Laplace’s equation over curves, surfaces, or regions by behav- 
ing like the potentials of distributions of charges over these re- 
spective loci, are similarly treated for reflection across x =0 by 
replacing the integrals by finite sums, reflecting the point charge 
singularity of each term, and passing to a limit. 


6. The Two-Dimensional Green’s Function Corresponding to 
(4). The two-dimensional case can be treated similarly to the 
three-dimensional one. The condition (8b) is replaced by the 
requirement that, near Po=(h, 0), u is to become infinite like 
—In ro= —In[(x—h)?+y?]"/2; the requirement (8c) at infinity 
might be kept, though the condition u = O(In r) in its place would 
be suggested by the theory of the logarithmic potential. One 
proves similarly that 


—h 
(23) u = —Inro—Inr, + 20 f in dz’ . 
The interpretation in terms of a positive image and an exponen- 
tial trail of negative images is obvious. 
For the present (two-dimensional case) the representation of 
u as the potential of a charge distribution is not at all unique. 
We shall prove this for real positive a. 
Indeed, by integration by parts and introducing imaginaries, 
(23) may be put in the form 


zt+h 


“ 


where z=x+7y* and “R” denotes “the real part of.” Replacin 
p g 


* No confusion is likely to arise from the use of z in a sense different from 
its previous one. 


= 
= 
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the (dummy) variable of integration x’ by z’, one may deform 
the path of integration from the real axis into an arbitrary curve 
C between the same limits. Thus 


sth eathtz") 
(24) “= R| in + 2f i], 
z—h c 2-2 
h dz’ 
+ 2f ( ) 
h c 


dz’ 
-2f ( 
2-2 


where J(z) refers to the “imaginary part of z.” Now the real as 
well as the imaginary part of 1/(z—2z’) represents a potential of 
a doublet or dipole placed at z=2’. Hence the curve integrals 
in (24’) can be interpreted as the potential of a proper distribu- 
tion of (logarithmic) dipoles over C, with axes respectively nor- 
mal to and along C. The tangential dipoles can by integration 
by parts be reduced to a distribution of poles. 

Of the various charge distributions the one derived from (23) 
is, no doubt, the simplest one. 

By integrating (24) one obtains 


s+ 
u= Rin 
(24’) 


sth 


where Ei denotes the “integral exponential.” The latter is multi- 
ple-valued. Another way of proving that u may be considered 
to be the potential of a variety of distributions is by rendering Ei 
single valued by using a branch corresponding to a cut along the 
curve C. The resulting discontinuities in C and its normal deriv- 
ative, when divided by 27, yield the densities of normal doublets 
and of poles respectively. 

The Fourier integral (in y) of u is also of some interest. Start- 
ing with the form (24) and utilizing 


1 
— = f for R(z)>0, 
0 


one obtains 
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+ h x 
in 2f a| 
z—h 


eh (ath) cos AY 
0 


(26) 


The Fourier integral as well as (25), (26) apply for R(@) <0, 
as well. 


7. Other Differential Equations. A pplications. As stated in the 
introduction, the reflection of singularities across a plane bound- 
ary corresponding to the boundary condition (3), in terms of a 
positive image and an exponential trail of negative images, can 
be carried over to certain other differential equatiozs. Among 
these are 


(27) — u = 0, = constant, 


(28) du/at = kV2u. 


The proof is again carried out by forming the function w given 
by (5), reflecting it, then solving for u. 

For (28), the heat conduction equation, the boundary condi- 
tion (4) represents a loss of heat (for real positive a) at the plane 
boundary x=0 at a rate proportional to the boundary tempera- 
ture. If we consider singularities due to instantaneous or perma- 
nent “heat sources” (while (4) applies), the temperature can in 
every case be determined by using a similar image source as 
well as an exponential trail of “sinks.” 

For steady point or line heat sources, (10) and (25) essentially 
still apply since (28) reduces to Laplace’s equation. As an ex- 
ample _of an instantaneous source recall that 


(29) v = for 0, = 0 for t <0, 


represents the temperature in an infinite medium due to an 
instantaneous line source of heat at the time ¢=0 along the 
z-axis, the amount of heat discharged being so chosen that 
| fedxdy=1 for t>0.* The temperature due to a similar line 
source at (h, 0) when (4) holds is given by 


* See Carslaw, The Mathematical Theory of Conduction of Heat in Solids, 
Chap. IX, 
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(30) 

The last integral can be expressed in terms of the “complemen- 
tary error” function defined by: erfc(x) 

Another application of the boundary condition (3) that will 
now be considered is in connection with eddy currents excited in 
a conducting semi-infinite solid x<0 by alternating currents 
flowing in free space in the region x>0. Time is assumed to 
enter as a factor e**‘. If angular frequency w is not too high, 
displacement currents may be neglected; then electromagnetic 
field theory, as expressed, say, by Maxwell’s equations, shows 
that the components of both the magnetic and electric field in 
the solid conductor satisfy the differential equation 


(31) ; V7u = au, 
where 
(32) a=(1+1)8, B=(2xwdy)!/”, \=conductivity, permeability. 


In free space the field components are harmonic in regions free 
from currents; near each current element di=(l, m, n)|i|ds, 
the electric and magnetic fields become infinite respectively like 


(33) E = 4riwdi/r, 
H = 4vX(di/r) 


a a 
tr|i|[(m= (+++), (Jam, 


where 7 is the distance to the current element, while /, m, n 
are its direction cosines. At the boundary x=0 the tangential 
field components H,, H.; E,, E, are continuous. 

We shall be concerned with the asymptotic behavior of the 
solution for large 8. It will be shown that for 6 large the differ- 
ential equation (31) can be replaced by proper boundary con- 
ditions on the components of the harmonic functions E, H. 

By examining cases where a complete solution is available it 
is found that for large B or large |a|, the terms 0?u/dy?, 0?u/dz? 
in (31) are negligible in comparison with the remaining terms, 
so that (31) approximately reduces to 


ll 


(34) 


= 
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(35) —_=a'n, 


whose pertinent solutions are 
(36) u = u(y, z)e%*. 


Physically, the constant 1/8 has the dimensions of a length, and 
when this length is small (compared to the dimensions of the 
exciting circuit and its distance from the plane x=0), the field 
developed in the conducting solid decays rapidly with penetra- 
tion into the solid and the resulting currents are confined to a 
thin layer near the boundary thus exhibiting the phenomena of 
“skin effect.” Using the above approximation we see then that 
within x<0 all the field components are of the form (36), 
where the coefficient of e** refers to their values at the boundary 
of the solid. Applying now the boundary conditions at 1 =0, one 
proves readily that the normal magnetic field H, satisfies the 
boundary condition 


0 a 
(37) =@Q at x=0. 
Ox 


Suppose that the current |i| flows along the closed curve Co. 
From (34) it will be seen that H,/4|i| becomes infinite at Co 
after the manner of the potential due to a distribution of 
doublets or dipoles along Cp and with axes parallel to the plane 
x=0. Applying the results of §5 to the reflection of H, across 
x=0 (the constant a is now equal to a/p and satisfies (9)), one 
is led (so far as H, is concerned) to an image of the original cur- 
rent along Ci, the reflection of Cy in x =0, as well as to a current 
sheet or solenoid along the cylinder which projects Co, C; on 
the plane x=0, the current density along the latter being 
—2(a/pye**=t" per unit x for x<—h. The remaining field 
components may be shown likewise to agree with the field of the 
original current and the images just described. 

The distributed currents, it will be noted, are not “in phase” 
with the original current. 

If Co is parallel to the plane x=0, then the electric field is 
likewise parallel to it, and the boundary condition (37) is satis- 
fied by its components E,, E,. 

It is possible to express the induced electromotive in terms of 
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the mutual inductances between Cy and its images in planes 
x =constant. It is intended to treat this subject at greater length 
at a future date and to present an exact treatment, including 
the nature of the reflections when (31) is not approximated by 
(35). 


8. Higher Normal Derivatives. In this concluding section we 
shall consider the extension obtained when (4) is replaced by a 
boundary condition involving higher derivatives: 


u 
(38) along x=0, 
x x 


where a and 6 are constants; a proper restriction on a and 6 will 
be indicated below. It will be shown that a singularity at P» 
corresponding to a point charge will be reflected into two ex- 
ponential trails along with a point charge. 

Perhaps the most elegant manner of arriving at this conclu- 
sion is by means of the algorithm of Heaviside’s operational 
calculus.* 

Denote the charge density of u along the x-axis by p(x) and 
replace 0/0x by p. At Po the charge density of u would be 
denoted by So(x—h), where the singular function So(x) repre- 
sents the “unit impulse” at x =0. Applying the left-hand opera- 
tor of (38) to u, one obtains a harmonic function which vanishes 
at x=0. The charge distribution of this function is given by 
(p?+ap+b)p(x), and for x>0 this reduces to (p?+ap+5) 
-So(x—h). The negative reflection in x=0 leads to the further 
charge — (p?—ap+5)So(x+h). Hence 


(p? + ap + b)p(x) = (p? + ap + b)So(x — h) 


3 
(39) — (p? + ap + b)So(x + A). 


Solving for p(x), we obtain 
*—apt+b 
p(x) = So(x — hk) 


So(x 


(40) 
A B 
= S§,(x — h) (1 + —— + +h). 
* The reader who is unfamiliar with the latter will prefer to regard the fol- 
lowing not as a proof but as a heuristic lead toward a result which is to be 


substantiated in some other way, for instance, by means of Fourier or Brom- 
wich integrals, or by direct substitution in (38) of the potential derived. 
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To obtain the last equation the partial fraction resolution 
ap+b 2a A B 

p p p ( ) 

+ap+b p> + ap+b 


was ued; here a, 8 are the roots of p?+ap/+5=0. Interpret- 
ing (40) for the case where R(a)>0, R(B)>0, one obtains 


p(x) = So(x — h) — So(x + h) 
-= + ]H(— x — h),* 


(41) 


(42) 


where //(x) is unity for positive x and zero for negative x. Thus 
there is a negative point charge at x= —h and a distributed 
charge of density —[Ae~7+%+Be“=*™] for x<—h. Since 
A8B+Ba=0, the total amount of distributed charge vanishes. 

It is of interest to point out that the content of (10), (11), 
and (22) can be derived in a similar fashion. In the former case 
one condenses the argument of §3 into 


(43) (p — a)p(x) = (p — a)So(x — h) + (p + a)So(x + h); 

in the latter case, the charge density at Po is pSo(x—h), and 
one obtains, by applying (p—a) and using negative reflection 
inx=0, 

(44) (p — a)p(x) = (p — a) pSo(x — h) — (p + a) pSo(x + hi). 
Solving (43), (44) for p(x) and using partial fractions, one ob- 
tains 

(45) p(x) = So(x — h) + So(x + hk) + 2a/(p — a)So(x + h), 
(46) p(x) = pSo(x — h) — [2a + p + 2a*/(p — a)|So(x + h), 
respectively, and, interpreting [1/(p—a)]So(x+A) as in (40), 
one is led to the charges implied in (10), (23). 


GENERAL ELEctrRic Co. 
ScHENEcCTADY, N.Y. 


* The present interpretation of [1/(p—a)]So(x+h) differs from that of the 
orthodox operational calculus and corresponds to using a Bromwich integral 
with the path of integration lying to the left of the origin. This could have been 
avoided by placing the original singularity along the negative x-axis. 

The more customary notation is / in place of H(x) and p/ in place of So(x). 
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THE CALCULUS OF VARIATIONS APPLIED TO 
NORLUND’S SUM* 


BY TOMLINSON FORT 


Néorlund’s sumf function 


Ss (x)Ax 


has many resemblances to a definite integral. The purpose of 
the present note is to point out how some of the classical meth- 
ods of the calculus of variations can be applied to such a sum. 
It may be that the field will prove fruitful for further research. 

We shall consider the problem of minimizing (maximizing) 
the sum 


where we have exactly the same understanding of what con- 
stitutes a minimum as in the classical problem of the definite 
integral. 


1. Euler’s Equation. We shall seek a necessary condition that 
a continuous real y minimize 


(2) < F(x, y, Ay)Ax 


similar to Euler’s equation for the corresponding integral. The 
condition of fixed end points in the integral problem is here re- 
placed by the condition that y be fixed over the interval 
Sc+1 and at the point 

For brevity in writing denote Ay by y’, F(x, y, Ay) by F(x), 
and assume that y is continuous and that F, F,, and F, are 
continuous in their arguments throughout all neighborhoods 
considered in the sequel. 


* Presented to the Society, September 9, 1937. 
+ Milne-Thompson, The Calculus of Finite Differences, page 201. In the 
present paper the difference interval is assumed to be 1. 


| 


886 TOMLINSON FORT {December, 


Assume that y minimizes (2). Give to y a variation en(x) 
where n(x)=0 when xSc+1, when x=b, and when x>M, 
where M equals } plus an integer. We assume also that n(x) 
is continuous throughout its range of definition. Let An=7n’. 

Under these restrictions 


b 
o(e) = S Fe, y ten, + e7’)Ax 
converges and is a function of e. 
We can differentiate under the integral sign inasmuch as all 


that is involved is an integral over a finite interval and a finite 
sum. We thus obtain the formula 


“(0 | - Sw. + 


nF + nF yAx. 


Summation by parts* yields 


b b z 
4 n (x)FyAx = POS n (x)Ax -f (x)Axdx 


c 


Fy-(b)n(6) — Fy (x)n(2)dx 
+ 1)AF,(x)Ax 


AG + 1)AF,.(x)Ax 


b+1 
n(x)AF,-(x — 1)Ax 


=— n(x)AF,-(x — 1)Ax since n(b) = 
e+1 


* Milne-Thompson, page 206. 
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Moreover, since n(x) =0 also when cSxSc+1, 


x)F ,(: Ax = x)F,(x)Ax. 
Consequently 

d b 

[F,(x) — AF, (x — 1)]Ax. 


This is necessarily zero for a minimum. Then necessarily 
(3) F,(x) — AF, («x — 1) = 0, x2ct+1. 


Suppose it were not true at #2c+1 where #+b+n, (n=0, 1, 
2,---). Assume 7=0 at all points except a neighborhood of 
Z not extending to any b+, and throughout which F,(x) 
—AF,(x—1) retains the same sign. For an 70 throughout 
this neighborhood this yields d¢/de |.-0~0. Moreover, since y 
is continuous it must also satisfy (3) when x=b+7; also M is 
as large as we please. Hence, a necessary condition that a con- 
tinuous y render (2) a minimum (maximum) is the satisfaction 


of (3), which we write 
(4) F,(x + 1) — AF, (x) = 0, 


The resemblance of this to Euler’s equation for the integral is 
immediate.* 
By exactly similar arguments making the assumptions that 


n(x) =0, cSxSc+n, and n(b)=n(b+1)= --- =n(b+n—1)=0, 


we arrive at the following analogue for Euler’s equation as a 
necessary condition for the rendering of (1) a minimum (maxi- 
mum): 


--- +(—1)"A*F =0, C, 
where y™=A"y. 


LEHIGH UNIVERSITY 


* Some special cases of the difference equation (4) were studied in detail 
by the author in the American Mathematical Monthly of March 1936. In 
that paper the independent variable was limited to integral values. Necessary 
modifications to make the discussions given there apply when the independent 
variable is continuous are not considered of sufficient interest to warrant 
printing. 
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ADDENDUM 
BY A. D. MICHAL AND E. W. PAXSON 


In our recent paper* on the application of mappings of 
topological neighborhoods on an open Banach subset 2 to 
abstract continuous group theory, we found Theorem 4 funda- 
mental. This was to the effect that every homeomorphic map- 
ping of a 2! on a (21, is an abstract coordinate 
transformation. We remark that a rigorous proof of this is as- 
sured if the topological space T is subjected further to the 
second Hausdorff countability axiom.f 

Instead of making use of this axiom we prefer to adopt ab 
initio an alternative mapping. With the notation of our paper 
we require the following interspace postulates to hold between 
T and the Banach space B. 


(1) For every Uc T, there exists at least one homeomorphism 
Uh, where = is an open set in B. Let K be the totality of such 
homeomorphisms. 

(2) All such sets © lie in a fixed open set 'c B. 

(3) There exists at least one Up ¢ T such that UAT. 

(4) If Uh isin K and U'c U, then U’h>’ isin K. 

(5) If isin K and then Uh>d’ is in K. 


Appealing to the continuity of Uhl, we may prove the 
following theorem: 


THEOREM 4’. Let Then there exist open sub- 
sets and ¢ D2) such that ts an abstract coordi- 
nate transformation. 

With this new mapping and with Theorem 4’ replacing 
Theorem 4, the formal deduction of the representative Lie 
equations proceeds as before with minor changes in interpreta- 
tion. 


CALIFORNIA INSTITUTE OF TECHNOLOGY 


* Maps of abstract topological spaces in Banach spaces, this Bulletin, vol. 42 
(1936), pp. 529-534. 

+ F. Hausdorff, Mengenlehre, 1927, p. 229. See also W. Sierpinski, General 
Topology, 1934, p. 40. 
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